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1 Introduction 

Amphiphilic molecules (the word derives from the Greek a/j,(j)l 4>l\Loi, meaning "love on both sides") are 
molecules which both love and hate water. They are formed by two parts with very different tastes, which 
are covalently bound together: one, the hydrophilic head, is polar or even ionized, and tends therefore to be 
close to the small, polar water molecules; the other, the hydrophobic tail, is usually a hydrocarbon chain, 
which perturbs the high order of water, and has therefore the tendency to pack close to similar chains. 

Amphiphilic molecules in solvents (water and/or oil) may form several structures (micelles, hexagonal 
phases, cubic phases. . . ), but we shall discuss mostly the cases in which they form a bilayer, i.e., a sheet 
made up of two layers of amphiphilic molecules: in water, the hydrophilic heads stem out of the bilayer on 
both sides, while the corresponding tails remain at the interior. Membranes can be formed by an isolated 
bilayer or by several bilayers stuck one on top of another: in this case one speaks of multilayers. 

Amphiphilic membranes are a physical realization of fluctuating surfaces. They are thin sheets (50-100A) 
of amphiphilic molecules immersed in a fluid, usually water or brine (water and salt). They can be of natural 
or artificial origin: the most important example of natural membranes is the cell membrane, which separates 
the interior of all living cells from its exterior. Living cells, and in particular eucaryotic ones, possess a 
large number of membranes, like the nuclear membrane, which separates the nucleus from the rest of the 
cell, allowing, e.g., mRNA to pass from the inside to the outside, and several chemical signals in the reverse 
direction, or the Golgi apparatus, which acts as a sort of "chemical factory" for the cell. Mitochondria are 
also organelles essentially formed by a membrane, folded on itself several times. 

Artificial amphiphilic membranes have recently become a lively research field stimulated by their appli- 
cations in the industry, in medicine and in cosmetics. One can form with them tunable or "active" filters, 
simulating, as it were, the action of the cell membrane. They are also able to close on themselves, forming 
vesicles (small closed surfaces), which may act as drug carriers, designed to open up and release their load 
when the "correct" physico-chemical conditions are found. Several other applications can also be envisaged. 

The physics of amphiphilic membranes is a wide subject, and it is out of question to review it fully in 
this series of lectures. I shall mostly dwell on the aspects which fit more closely the scope of the School, 
namely those involving shape fluctuations. It will be necessary to review briefly the basic physical chemistry 
involved to understand why membranes form at all, which features govern their shape, and their equilibrium 
or dynamical behavior. General introductions to the statistical mechanics of amphiphilic membranes can 



be found in ref. |93|]. A goo d introduction to the basic properties of biological membranes is found in [114 



Chap. 12]. References 109 , pi| contain introductions to the physical chemistry of amphiphilic molecules. 

Therefore, in the next section I shall briefly dwell on the structure of their basic components, the am- 
phiphilic molecules, and give an overview of the structures they form in the presence of water and/or oil. In 
the following section I shall discuss the free energy of an isolated fluid membrane as a function of its shape. 
The corresponding hamiltonian (due to W. Helfrich) lies at the basis of the current understanding of vesicle 
shapes and of membrane fluctuations. The following section contains a brief review of recent theoretical 
and experimental work on vesicle shapes. Then the effects of fluctuations on amphiphilic membranes will be 
discussed: these involve on the one hand the characteristic flicker phenomenon in vesicles, and on the other 
hand the renormalization of the elastic parameters appearing in the Helfrich hamiltonian. All these aspects 
are reflected in the phase behavior of interacting fluid membranes to which the last section is dedicated. A 
few technical points are discussed in the four Appendices. 



2 Amphiphilic molecules and the phases they form 

Biological membranes are formed by a bilayer of amphiphilic molecules, the most common of which are 
phospholipids. Their hydrophilic head is a phosphate, and their tail is formed by one or two fatty acids. Most 
often the two parts are connected by a glycerine hinge, and the molecule is therefore called a phosphoglycerid. 
The chemical structure of a typical phosphoglycerid is represented in fig. |]. 

The glycerine molecule which forms the hinge of the structure is shown on the top. On one side, it is 
ester- linked to a phosphate group -POOH-O-R carrying the ethanolamine residue R=-CH2-CH2-NH2. 
In physiological conditions the head is almost always ionized, yielding -PO _ -CH2-CH2-NH^. Via the 
other two carbons (1 and 2) the glycerine is ester-linked to two fatty acid chains, of the general struc- 
ture CH3-(CH2)„-COOH. In our case one has n = 10, and the acid is called acid n-dodecanoic, or lauric 
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Figure 1: Chemical structure of a typical phosphoglycerid: 1,2 dilauryl-DL-phosphatidylcthanolaminc 
(DLPE). From B. L. Silver, The physical chemistry of membranes (Boston: Allen & Unwin, 1988) p. 2. 

acid. The example helps to clarify the terminology. Other phosphoglycerids may differ from DLPE either 
by having just one fatty acid chain (monoglycerids) , or by the nature of the fatty acid tail. For example, 
if the two chains have n = 14 (palmitic acid) one has dipalmitoylphosphatidylethanolamine, (DPPE). If we 
change now the residue R to choline Ri =-CH2-CH2-N + -(CH3)3, we have dipalmitoylphosphatidylcholine 
(DPPC), also called lecitin, which is an important component of natural membranes. 

In general, the two fatty acids of natural phosphoglycerids are different, and one of them is quite often 
unsaturated. The knee which appears because of the insaturation, like that shown in fig. ||, helps in keeping 
the fluidity of the membrane. 

Pure phospholipids may be made to crystallize, and the crystal structure of several of them have been 
determined by X-ray or electron diffraction. In the crystalline state they are stuck in bilayers, their tails all 
trans, with their heads folded approximately parallel to the bilayer surface and linked (mostly by hydrogen 
bonds) into a firm network. This organization is schematically represented in fig. ||. This order is disrupted, 
as the temperature increases, via a two-step process: 

• At a temperature T t called the transition temperature the order of the chains breaks down: i.e., a 
significant fraction of the carbon links goes over to the gauche configurations (cf. fig ||) , providing a 
gain in entropy against a loss in van der Waals attraction among the chains. This is sometimes called 
"chain melting" or "premelting" . 

• The actual melting temperature T m is reached when the ionic lattice formed by the hydrophilic heads 
eventually breaks down. 

The melting temperature T m of most pure phosphoglycerids is comparatively high (~ 200°C), and does not 
depend strongly on the length of the fatty acid tail. On the other hand the transition temperature is closer 
to room temperature and increases with the hydrocarbon chain length: a plot is shown in fig. |5[ We expect 
therefore, in most phospholipids, to find a phase in which a relatively fluid hydrocarbon core is sandwiched 
between two relatively rigid polar sheets. This feature makes possible to sustain fluid bilayers in an aqueous 
medium. In isolated bilayers in a solvent one observes a sharp anomaly of the specific heat at a transition 
temperature T t somewhat lower than that observed for pure anhydrous phospholipids (see fig. ^). This 
anomaly is accompanied by a fast variation of the mechanical properties of the bilayer, and is known as the 
"main transition" . 

The idea to keep in mind is that amphiphilic molecules are complex ones, with a number of internal 
degrees of freedom in the tails. The transition takes place essentially in their tails. The disorder in the tails 
eventually induces a liquid-like disorder in the location of the molecules on the sheet. This is quite different 
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Figure 2: Space filling models of (A) palmitate (Ci6, saturated) and (B) oleate (Cig, unsaturated). The cis 
double bond in oleate produces a bend in the hydrocarbon chain. From L. Stryer, Biochemistry (New York: 
Freeman, 1988) p. 285. 



from a melting transition in a sheet of bead-like molecules. The coupling between crystallization and chain 
ordering has not yet been satisfactorily treated |5| . Its relevance can be grasped from the observation that 
most biological membranes appear to work very close, but slightly above the transition temperature. Some 
suggested explanations of this fact can be found in H . 

As soon as water enters the picture, it makes it much more complicated. With a low water content, 
the system maintains by and large the lamellar organization characteristic of the pure crystal. The layers 
can exhibit a number of different phases, which have different stability domains. For example in a mixture 
of dimyristoylphosphatidylcholine (DMPC) and water, one can find at least three phases at high DMPC 
concentration |6^, 4(| |^J. They are called the L Q , Lp and phases. In the L a phase the bilayers arc 
fluid and flat on average. If one lowers the temperature T, or decreases the water content <j>, one goes to an 
ordered, "solid like" phase L^, in which the hydrocarbon chains are ordered and the molecules do not diffuse 
freely. In these phases, the order of the hydrocarbon chains implies a larger thickness of the bilayers. One 
can also observe an intermediate "rippled" phase Pg (see fig. 0), in which the bilayers exhibit an undulated 
structure and almost solid- like diffusion properties fl67[. Analysis of X-ray experiments on these "rippled" 



phases strongly suggests that they are characterized by a modulation of the bilayer thickness |118|. The 
hydrocarbon chains are often tilted with respect to the bilayer: one then denotes the phases as ~Lp> or Ppt. 
In fact, there are several different L^/ phases, distinguished by the relationship between the tilt and in-plane 
bond orientational order [112|. 

The "rippled" P^ phases may be considered as an intermediary structure between the lamellar ones found 
at low water concentration and the ones found at low amphiphile concentration. If we add amphiphilic 
molecules to pure water, the molecules first go preferentially to the air- water interface, forming a monolayer, 
with their heads toward water as long as the concentration does not exceed the critical micelle concentration 
(cmc), which is of the order of 10 _10 mol. Below the cmc the amphiphilic molecules are overwhelmingly in 
the monomer form, at higher concentration added monomers appear almost exclusively in aggregates, mainly 
of globular form with the hydrophilic heads on the surface. These aggregates are called micelles. They form 
more readily for single-chain amphiphiles (e.g., monoglycerids) and are favored by the presence of large head 
groups. The structure of micelles is depicted in fig. ^. In this figure, the concentration of amphiphiles is 
high enough to let the micelles arrange in a close-packed bec lattice: this is the simplest example of the 
remarkable cubic phases formed by amphiphiles. As the amphiphile concentration is increased, one observes 
the appearance of nonspherical micelles, and eventually of cylindrical rods. These rods behave as "living 
polymers" at low concentrations, and organize as a close packed, hexagonal phase, like that shown in fig. ^, 
at higher concentration. 
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Figure 3: Crystal structure of DMPC projected onto the crystallographic a-c plane. 
I. Paschen, Nature 281 499 (1979). 



From R. H. Pearson, 
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Figure 4: Conformation of C-C bonds in fatty acyl chains: (A) trans (t) conformation; (B and C) a 120° 
rotation yields a gauche (g) conformation, which can be either g + (clockwise rotation) or g _ (counterclockwise 
rotation). From L. Stryer, Biochemistry (New York: Freeman, 1988) p. 297. 
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Figure 5: Chain melting temperatures T t for some amphiphiles as a function of chain length. + anhydrous 1,2- 
diacyl-DL-phosphatidylethanolamincs; •, 1,2-diacyl-DL-phosphatidylethanolamines in water; A, anhydrous 
1,2-diacyl-L-phosphatidylcholinc monohydrates; o, 1,2-diacyl-L-phosphatidylcholine monohydrates; x, 1,2- 
diacyl-L-phosphatidylcholinc monohydrates in water. For comparison, the melting points of some normal 
paramnes, □, are given. From D. Chapman, et al., Chem. Phys. Lipids 1 445 (1967). 
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Figure 6: Differential scanning calorimetry of phosphatidyl choline. The small peak (the pretransition) 
comes from a change in the tilt of fatty acyl chains with respect to the bilayer. The major peak arises from 
the phase transition in which crystalline fatty acyl chains become disordered because of the introduction of 
kinks. From L. Stryer, Biochemistry 3rd ed. (New York: Freeman, 1988) p. 297. 




Figure 7: Electron microscope photograph of "rippled" phases in DPPC, produced by cooling large vesicles 
from the P^' phase. From E. Sackmann, et al., Liquid crystals of one- and two-dimensional order, (Berlin: 
Springer, 1980) p. 314. 



Figure 8: Scheme of micellar structure. The spheres indicate the head groups and the wiggly lines the lipid 
chains of phospholipids. From C. L. Khetrapal, et al., Lyotropic liquid crystals (Berlin: Springer, 1975). 



If we add paraffinc oil, close in composition to the hydrocarbon tails of phospholipids, we may stabilize 
some new phases. A schematic picture of the resulting phase diagram is found in ref. |}2| (fig. [To]). The new 
phenomena are determined by the fact that, if the oil content is large enough, amphiphilic molecules can 
form a monolayer, with their tails towards the oil-rich phase and their heads towards the water-rich one. 

This opens the way to bicountinuous phases, in which both the oil-rich and the water-rich phases per- 
colate. Some of these phases have a periodic, cubic structure, similar to that represented by the cartoon in 
fig. the actual structures are however more complex and difficult to draw. They are collectively known 
as "plumber's nightmare" phases, since they may be considered as pipeworks in which the interior and the 
exterior look the same. Of great interest, both theoretical and experimental, are the irregular bicontinuous 
phases. As shown from the phase diagram, they can be obtained by adding oil to a water-amphiphile solution 
(with a high enough amphiphile content), without crossing any phase barrier. Moreover they can coexist 
with the water-rich and the oil-rich phases, having a water-oil ratio close to one. In this coexistence they will 
remain between the water and oil phases because of their density, and are known therefore as middle-phase 
microemulsions. The name microemulsions intimates that they are formed of almost equal proportions of oil 
and water, like ordinary emulsions: however, whereas emulsions are nonequilibrium structures obtained by 
suspending droplets, say, of oil in water by means of intense mixing, microemulsions are equilibrium phases. 
In nonequilibrium emulsions, the droplet size ranges from the micrometer to a fraction of millimeter, whereas 
in microemulsions one observes irregularities in the local composition at scales of the order of a few hundred 
Angstroms, much smaller than the wavelength of visible light. As a consequence, microemulsions usually 
appear transparent. 

There is also the rather paradoxical possibility that the bilayer forms an interface separating two perco- 
lating domains occupied by the same solvent. This is known as the sponge phase [ 105 1 . 

The structure of these phases can be studied by the "traditional" means of light or X-ray diffraction 
(depending on the characteristic size of the structures), but can be also directly exhibited by freeze fracture. 
In this technique the sample is rapidly frozen to the temperature of liquid nitrogen. The frozen sample 
is then fractured by means of a microtome knife. Cleavage usually occurs in the middle of bilayers. The 
exposed regions can then be shadowed with carbon or platinum, which produces a replica of the interior of 
the bilayer. In order to expose the exterior of the membrane, one can combine freeze fracture with etching. 
First, the interior of a frozen membrane is exposed by fracturing; then, the ice that covers one of the adjacent 
membrane surfaces is sublimed away: this process is called deep-etching. The combined technique, called 
freeze- etching electron microscopy, provides a view of the interior of a membrane and of both its surfaces. 
The structure of emulsion and microemulsions can also be exhibited in this way. 

By use of these techniques it has been possible to prove the general validity of the fluid mosaic model of 
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Figure 10: Schematic oil-water-surfactant phase diagram with microstructures depicted. From H. T. Davis 
et al., in J. Meunier, D. Langevin, N. Boccara (eds.), Physics of amphiphilic layers (Berlin: Springer, 1987) 
p. 311. 
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Figure 11: Freeze-fracture of a bicontinuous microemulsion containing equal volumes of water and oil. The 
shadow material decorates the oil-rich parts of the fracture face. Bar: 500 nm. From W. Jahn, R. Strey, in 
J. Meunier, D. Langevin, N. Boccara (eds.), Physics of amphiphilic layers, (Berlin: Springer, 1987) p. 354. 




Figure 12: Freeze-fracture of a cell membrane from an onion root tip. The uppermost area in the photograph 
is the cell cytoplasm with a vacuole (V) evident. The center area shows the fracture face of the cell membrane 
(plasmalemma, PL). The fracture has also ruptured the cell wall (CW). Plasmodesmata (P), fine threads 
connecting the cytoplasm of adjacent cells, have also been broken off. From D. Branton, D. W. Deamer, 
Membrane structure (Berlin: Springer, 1972). 
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Figure 13: Schematic illustration of a biological membrane (courtesy of Ove Broo S0rensen of the Danish 
Technical University). The membrane itself is a mixture of several kinds of phospholipids, plus smaller am- 
phihilic molecules, among which (in animals) cholesterol has a prominent role. The dark objects traversing 
the membranes are "integral membrane proteins" . Other, more flexible, membrane proteins carry carbohy- 
drates (represented as small polygons) in their exterior. The rubber-like cytoskeleton attached to the inner 
side of the membrane exerts a strong influence on the mechanical properties of the composite membrane. 



cell membranes, proposed by S. J. Singer and G. Nicholson in 1972 [111 . It may be schematically represented 



by the cartoon in fig. [12|. The main part of the membrane is formed by a bilayer, which is a mixture of 
several kinds of amphiphilic molecules, in particular phospholipids and glycolipids (sugar-containing lipids, 
like in particular sphingomyeline and cerebroside), plus smaller ones like cholesterol. Integral membrane 
proteins are dissolved in the bilayer: they can freely diffuse laterally, but cannot move out of the surface. 
Other proteins carry hydrophilic tails (sometimes containing sugar: one then speaks of glycoproteins) which 
extend in the exterior of the cell. In the interior of the cell there is often a network of filaments, suitably 
anchored to the bilayer. This is the case in particular of the red blood cells (erythrocytes), whose skeleton is 
formed by filaments of spectrin bound to "buoys" formed by other proteins, like ankyrin and "protein 4.1". 

The fluidity of biological membranes can be exhibited by fluorescence photobleaching recovery experi- 
ments in intact cells. One first attaches a fluorescent dye to a specific membrane component. One then looks 
at a small region (~ 3/im 2 ) through a fluorescence microscope. One then destroys the fluorescent molecules 
in this region with a very intense light pulse from a laser. The fluorescence of this region is then monitored 
as a function of time. The rate of recovery is related to the diffusion coefficient D of the fluorescent-labeled 
molecules. The order of magnitude of D turns out to be 10 -8 cm s _1 , implying that a molecule can diffuse 
about 2/im in 1 s. On the other hand, the characteristic times for a molecule to pass from one membrane 
layer to the other (flip-flop) are of the order of several hours. These times can be measured, although with 
some difficulty, by NMR J74|. Therefore any asymmetry of the lipid bilayer can be preserved for long periods. 

Of course a biological membrane is extremely complex to characterize physically. A strong tendency has 
therefore developed towards the study of model membranes, containing just one (or two) phospholipids, and 
in case a controlled concentration of impurities. For example, one can suspend the amphiphile in an aqueous 
medium, and the agitate the mixture by high-frequency sound waves. This procedure is called sonication. 
Alternatively, on can dissolve the lipid in ethanol, and then inject the solution via a fine needle into water. 
In this way one obtains aqueous compartments closed by a lipid bi- or (more often) multi-layer. These 
structures are called lipid vesicles or liposomes. It is of course possible to prepare the liposomes in a solution 
containing some interesting drug, and then separate them from the surrounding solution by dialysis or by 
gel filtration. This technique provides in principle a way to control the delivery of drugs to target cells. 
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3 Isolated membranes: the Helfrich hamiltonian 

In this section I shall discuss the fundamentals of the physical description of isolated membranes: I shall 
neglect therefore the interaction between membranes, and I shall not consider the microscopic mechanism 
leading to their formation. 

We can write the free energy of a bilayer composed of 2N amphiphilic molecules immersed in water as a 
function of the area a — A/N per molecule 

F = 2N(f>{a) = 2N [0 phob (a) + ph n(a) + &nt(a)] . (3.1) 

In this equation, <p{a) is the chemical potential of molecules in the bilayer, p hob and p hii are the effective 
attractive and repulsive parts, respectively, due to their interaction with water, while 0j nt is a direct interac- 
tion between amphiphiles which will generically be repulsive. In general, 4>(a) will have a minimum at some 
preferred area-per- molecule a = oq. If the membrane does not exchange molecules with the reservoir, and if 
it can freely adjust its total area, the equilibrium will be reached when 



OF 
da 



-0, (3.2) 



implying 

a cq = a . (3.3) 
Therefore, the surface tension 7(a) of the membrane vanishes at equilibrium: 

OF 



7(« cq ) dA 



da 



0. (3.4) 



This simple thermodynamic argument is usually evoked to explain why the surface tension of amphiphilic 
membranes can be very small. 

If we consider a fluctuating amphiphilic membrane, it is necessary to distinguish between its total area 
A and its projected area A p |3jJ |. To fix one's ideas, suppose that the membrane spans a planar frame of 
area A p . If the compressibility of the amphiphilic monolayers is low, the total area A is proportional to the 
total number 2N of molecules forming the membrane: A = aoN. It can only vary through changes in the 
number of molecules N. Then the two quantities A and A p can be considered as independent thermody- 
namic variables fl(i|| . They are both extensive: their thermodynamic conjugates represent distinct physical 
quantities. The area coefficient conjugate to the total area A, which we denote by 7, is for incompressible 
fluids directly proportional to the chemical potential [i of amphiphilic molecules. The film, tension conjugate 
to the projected area A p , which we denote by r, corresponds to the physical "surface tension" . One can then 
consider four different thermodynamical ensembles: 

1. (A, j4 p )-ensemble: isolated, framed membranes. 

2. (A, r)-ensemble: isolated, unframed membranes. 

3. (7, ^4 p )-ensemble: open, framed membranes. 

4. (7, r)-ensemble: open, unframed membranes. 

Experimentally, the most important situations are the open, framed systems (which can be experimentally 
realized in black lipid films or monolayers) and the isolated, unframed systems, in which the projected area 
can fluctuate, which correspond to the case of lipid vesicles (at least as long as exchange of lipids with the 
surrounding solution can be neglected). 

Open, framed systems. 

In this ensemble the projected area A p is fixed and the total area A fluctuates. The fluctuations are governed 
by a hamiltonian of the form 

H = jA + H ch (3.5) 
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where H. e \ contains the contribution of elastic internal forces (bending energy, shear modulus, etc.). The 
partition function is written as the sum over all film configurations C with fixed A p : 



(3.6) 



c 

The free energy in this ensemble is 

Go (7,40 = ~kBT In Z ( 7 , A p ), (3.7) 
and the film tension is simply defined as the free energy per unit projected area: 

T= lim ^Al. (3.8) 

Ap^oo Ap 

Isolated, unframed systems 

In this ensemble the total area A is fixed while the projected area may fluctuate. The thermodynamic 
potential is obtained from Go by first going to the isolated, framed (A, A p ) ensemble where both A and A p 
are fixed, and the thermodynamic potential is obtained by a Legendre transform 



F I (A p ,A) = Go(A p ,r)-~,A, A ' " 



d~; 



(3.9) 



Then one goes to the isolated, unframed film ensemble by a second Legendre transform which defines the 
associated thermodynamic potential 

Gi(t,A) = F i (A p ,A)-tA p , (3.10) 
where the surface tension r is defined by 



dFi 



dA p 



dGr 



dA p 



(3.11) 



It is clear that in the thermodynamic limit, A — > oo, the surface tension defined in this ensemble by eq. (3.11) 
coincides with the surface tension defined for the open, framed system by eq. ( |3.9| ). 

We can now analyse the meaning of the tension r for isolated, unframed systems with fixed total area 
A. Let us assume that the projected area fluctuates around its mean value (A p ) . This mean value can be 
obtained by minimizing Fj(A p ,A) with respect to A pi while A is fixed. In the thermodynamic limit it is 
useful to consider the free energy density 

/ = (3-12) 

as a function of the area ratio 

a p = ^. (3.13) 

Two situations are then possible, as depicted in fig. |lj. 

• /(dp) has its minimum for a nonzero ratio < a p < 1 (fig. |l4|(a)). The membrane is then said to 
be flat, since, although shrunk due to thermal fluctuations, it still keeps the global structure of a 
two-dimensional object. In this case, eq. (|3.11 ) implies that the surface tension vanishes. 



/(a p ) has its minimum at a p — (fig. [l4|(b)). The membrane is then said to be crumpled, since it is 
so shrunk by thermal fluctuations that its extension in space does not scale linearly with its internal 
tension. In this case it is obvious from fig. |l4|(b) that the tension r has in general no reason to vanish: 
however, there might be marginal situations in which, although the minimum of f(a p ) is at a p = 0, its 
slope at this point vanishes, and one has therefore r = 0. 
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Figure 14: The free energy density / as a function of the area ratio a p . (a) / has a minimum for a p > and 
the tension r vanishes, (b) / has its minimum at a p = and the tension r may either be positive (i) or zero 
(ii). 



Let us consider a fluid membrane freely fluctuating in a solvent. We assume that its thickness is much 
smaller than the length scales describing its shape and its undulations. Its fluidity implies that all the internal 
degrees of freedom, related, e.g., to the hydrocarbon chain conformation, the local molecular density, etc., 
reach equilibrium on a fast time scale. We are thus led to the conclusion that the free energy of the membrane 
depends on its shape alone. In our hypotheses, we may describe its shape as a geometrical surface. 

The curvature of a surface can be described by two quantities, the mean curvature H, and the Gaussian 
curvature K . Consider a portion of the surface around an arbitrary point O, like the one depicted in fig. 
Draw the tangent plane through O and choose the local coordinate axes (t 1 , t 2 ,t 3 ) so that O is their common 
origin, t 1 and t 2 lie on the plane, and t 3 = h is normal to it. The surface is locally represented by a quadratic 
form: 

i 2 

h = -^2 n^v + ... (3.i4) 

The matrix f2 = (Oy) obviously vanishes for a locally flat surface. This expression of the distance h from 
the tangent plane in terms of Cartesian coordinates on it, is also called the second fundamental form of the 
surface. The curvature is defined via the invariants of the second fundamental form, namely: 

• The mean curvature H is given by H = TrO; 

• The Gaussian curvature K is given by K = DetSl. 

The eigenvalues c%, C2 of Q are called the principal curvatures. They are the inverse of the curvature 
radii of the intersections of the surface with two planes containing the normal to the surface, and mutually 
perpendicular. The intersections of these planes with the tangent plane define the principal curvature 
directions. The principal curvatures are the extreme values of the curvature of any intersection of the 
surface with a normal plane. In the exceptional case where c\ = c%, the point O is called an umbilical point, 
and the principal directions are not defined. The mean curvature H — c\ + c-i is positive (negative) if the 
surface is locally mostly above (below) the tangent plane, going in the positive direction along the normal. 
The Gaussian curvature K = c\pi is positive if the surface is locally on one side of the tangent plane (elliptic 
point), and negative if it is on both sides (hyperbolic point). 

The curvatures H and K have a simple geometrical interpretation. Consider a small portion S of the 
surface, and let AA be its area. Construct now a new surface S' by displacing each point of S of a distance 
S along the normal, in the direction of positive h. The new surface will have the area 

AA' = AA (l + SH + S 2 K + o(<5 2 )) . (3.15) 

We now make the following assumptions on the free energy of membranes: 

• The surface is smooth, and can be locally represented by the parametric equation r = r(a), where 
r denotes a point in three-dimensional ambient space, and a = (a 1 , a 2 ) are local coordinates of the 
surface. The function r(er) is differentiable an arbitrary number of times. 
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Figure 15: The tangent plane at a point O. The distance of a generic point P from it is given by h = 
I Tlij ^ijt 1 ^ ■ The principal directions on the surface at the point O are also drawn. 

• The free energy can be expressed as a local functional of r(a) and its derivatives. This assumption 
rules out, for the time being, the effects of the interaction of the membrane with itself, due for example 
to close contacts when a part of the membrane folds on the rest. 

• The free energy must be invariant under Euclidean transformations applied to r and under reparametriza- 
tion transformations like a — * a' = &'(a). 

It was noticed by Canham jlij and then by Helfrich |55| that these hypotheses imply that, considering 
only the contributions of derivatives of f up to second order, one obtains the following general expression of 
the free energy of a membrane: 



F= I dA 

s 



l+\n{H-H f + kK 



(3.16) 



Here dA is the area element and 7 is the area coefficient. The integral is extended over the membrane surface 
S. The coefficients k and k are known as the rigidity and Gaussian rigidity respectively. The quantity Hq 
is called the spontaneous curvature. 

The free energy fl3.16| ) has been derived only on the basis of geometrical considerations. However it is 
necessary to gain some insight on the terms appearing in this expression by considering in more detail the 



elasticity of a fluid membrane. One finds several discussions of this subject in the literature |99, 116, |57[,|8C|, 
A PP . B]. 

The simplest derivation starts by considering a monolayer. Let us assume that the elastic energy per 
molecule in a given configuration is the sum of the energy £^ pertaining to the heads and the energy £ t 
pertaining to the tails. We take as a reference the neutral surface, where the moments of the elastic forces 
on heads and tails cancel. We assume for definiteness that the positive direction is towards the heads, and 
denote by 6^ and 5t the distance of the heads and the tails, respectively, from the reference surface. 

The areas of the heads and at of the tails are given by: 

a h = a(l + H5 h + K5l + ...), (3.17) 
a t = a (l - H5 t + K5% + ...), (3.18) 

where a is the area of the molecule on the neutral surface (cf. fig. |l6| ). We can now write, assuming a simple 
(Hookc-like) elasticity: 

ft, = M^^y, (3.19) 



2 V aho 

ft = \h( a -^)\ (3.20) 

where aho and ato are the equilibrium values of head and tail areas, and and k% are elastic constants (with 
the dimensions of an energy per molecule). The elastic energy £ per unit molecule can be written 

£ = £ h +£t- (3.21) 
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Figure 16: Simple model of a bilayer in an aqueous medium. One assumes that the stresses are localized 
on the heads and tails of the amphiphilic molecules. The position of the neutral surfaces is denoted by the 
broken lines. 

For a flat membrane we can obtain from it the expression of the elastic energy of compression per unit area: 

Htlv = - = ^ + \k^( a -^)\ (3.22) 
a Q a Q 2 \ a J 

where Eq is a constant, and where the compression modulus fc*- m ) and the equilibrium area per molecule ao 
are respectively given by 

fc( m ) = ^L + A (3.23) 



«h0 a t0 



0, 



, a h0 a t0 / Ph + Pt 



The index (m) reminds us that we are dealing with a monolayer. The condition that the moments of the 
compression forces vanish on the neutral surface implies OhAi = o~tPt- We can now let the expression of the 
areas a t and ah in the expression of the elastic energy, obtaining 



■r-(m) _ oj(m) , ^/(ro) 
J ' "-comp i ' L bcnd 



(3.25) 



where 



a (m) = — , (3.26) 
a 

^ (m) = (S h + 5 t ) 2 k^—^—=S h S t k^\ (3.27) 

(Ph + Pt) 

_(m) _ / , PhPt(Ph~Pt) ,„ „„s 

" Ko-a t o) G g t + /3t)a , (3-28) 

= - r ^ T - (3.29) 

We see that K< m ) > 0, while 7t (m) can be of either sign. We can now estimate the order of magnitude of the 
rigidities For pure phospholipids the measured compressibility is k ~ 50 dyn/cm = 510 1 Jm , 

and the distances o"h, o"t are of the order of 10 A. We obtain therefore ~ 5 10~ 20 J, which must be 

compared with the termal energy k-gT ~ 4.210~ 21 J at T = 300 K. The fact that typical rigidities are 
only slightly larger than the thermal energy makes fluctuations important in understanding the behavior of 
amphiphilic membranes. 

We can now consider a bilayer, by putting two monolayers, called "internal" (i) and "external" (e) on 
the top of each other. We assume that there are no interactions between the monolayers, so that 

T h =^ (m) +^i m) , (3.30) 
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where the suffixes refer to the two monolayers. We denote by dj the area per molecule of the monolayer 
i = i, e, measured on its neutral surface. We can thus define the particle density pi of each monolayer, 
measured along a bilayer reference surface, placed midway between the two neutral surfaces of the monolayers. 
Let us denote by 6 the distance between the bilayer reference surface and that of the monolayers. For a 
bilayer with heads outside, for example, we have S = St- We choose the positive direction on the normal, in 
going, e.g., from the internal to the external layer. We then have 

' ai (l- H5 + K5 2 + ...), (3.31) 



Pi 
1 



= a e (l + H8 + K8 2 + ...). (3.32) 



Po 

This allows us to define the mean density p and the density difference p: 

P = (3-33) 
9 = (3-34) 

In the case where the spontaneous curvature H of each monolayer vanishes, one obtains the following 
expression of the elastic energy density per unit area of the bilayer: 

?h = cTbP«o + Tjfcb (aoP - I) 2 + ^k h ( a op) 2 
( h 5 2 \ 

+a k h DpH + ( 2K b + -y- J H 2 + K h K. (3.35) 

In this equation, the bilayer parameters (denoted by the suffix (b)) are simply the sum of the corresponding 
parameters of the monolayers. 

If the spontaneous curvature of the monolayers does not vanish, but the two monolayers are made of 
the same material, one obtains the same expression, but with slightly rcnormalized values of the elastic 
parameters. In either case, the spontaneous curvature Hq of the bilayer vanishes. 

The only difference with respect to the Helfrich hamiltonian lies in the terms depending on the particle 
densities p, p, which we now discuss. The first term, erbpao, yields a contribution proportional to the total 
number 2N of molecules in the membrane, but independent of its geometry. It corresponds to a shift 
in the chemical potential of the amphiphilc. As long as one considers vesicles over short time scales, so 
that amphiphilic molecules arc not exchanged with the solvent, this term plays no role. The second term 
represents the energy cost necessary to induce local density fluctuations. As we have seen, the rigidity Kb is 
of the order of fcb<5iA: therefore, it costs about the same energy to bend the membrane with a radius of the 
order 10 3 <5 <~ 1 /im as to produce a change of order 10~ 3 in the relative area (A — Aq)/Aq, where Ao = Nao. 
From a practical point of view, therefore, the membrane can be considered incompressible and its area can 
be set equal to A . The third term represents a coupling between the local density difference between the 
two monolayers and the mean curvature H. Indeed, if we try to bend the membrane, the molecules contained 
in the interior layer are compressed with respect to those contained in the exterior layer, and therefore tend 
to escape. 

Let us now consider a closed vesicle, formed by a bilayer with identical monolayers. If the observation 
times are not too long (of the order of several hours, or even a few days) the total number of molecules 
contained in the bilayer remains constant: on the other hand, as we have seen, the "flip-flop" times are also 
quite long, and therefore the number of molecules contained in each monolayer can be considered constant. 
Therefore the integrals of both p and p over the whole bilayer surface are constant. In particular, that means 
that their average values are given by 

Po = — (3-36) 

Po = — 2^4— i (3-37) 

where N c ^ are the number of molecules in the two monolayers, and A is the area of the vesicle. Therefore 
the term § dA-ypao is independent of the bilayer shape. It is reasonable to assume, in view of the fast lateral 
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diffusion times of the amphiphiles, that p locally equilibrates faster than the vesicle shape. It will then 
assume locally the value dictated by the optimization of Fy, when the curvatures are kept fixed. One has 

p(f)=p + Sp(f), (3.38) 

where Sp satisfies 

Sp = -po H. (3.39) 

a 



We can use again eq. ( |3.15| ) to show that this equation implies 

dASp = 0. (3.40) 



i 



Js 

If we assume that ~p — Oq 1 locally, and let the result for dp into the expression for Tt, we obtain 

d.U-i, = f <L i <j ' 



F b = j> dAF b = $ dA \ ^n h H 2 +K h K\, (3.41) 



up to a term independent of the bilayer geometry. 

These arguments lead thus to the conclusion that vesicles formed by symmetric bilayers are described 
by the Helfrich hamiltonian with a vanishing spontaneous curvature. A small asymmetry between the 
monolayers can be induced by a difference in the solvent: e.g., if a higher concentration in ions screens 
more effectively on the interior the electrostatic interaction among the heads. One is safe to assume that 
any induced spontaneous curvature will be rather small. However, the surface must satisfy the additional 
constraints, fixing the value of (i) the enclosed volume V, (ii) the total area A, and (iii) the area difference 
between the two monolayers, which can be expressed in terms of the integral of the mean curvature H: 

AA = 28 j> dAH = 25M. (3.42) 



This model of a vesicle, due to Svetina et al. [115], is known as the bilayer coupling (BC) model. The 



earlier work on vesicle shapes had considered the spontaneous curvature as an independent parameter, 
and had neglected the constraint on the integral of the mean curvature. This model has is known as 
the spontaneous curvature (SC) model. Although the spontaneous curvature Hq can be considered as the 
conjugate variable to the integral of the mean curvature M, the two descriptions are not equivalent. The 
point to keep in mind is that the Helfrich curvature energy does not grow with the size of the system as an 
extensive free energy, in the way one is accustomed to find in thcrmodynamical systems. In fact, it is scale 
invariant: consider a surface S defined by some parametric equation of the form r = rb((r). Now define the 
new surface S' by the parametric equation f = Xro(a), where A > 0. In this transformation, we have 

H —> H' = X^H- (3.43) 
K -► K' = \- 2 K; (3.44) 
dA -» dA' = X 2 dA. (3.45) 

Therefore the elastic energy remains locally invariant: on the other hand, since H (and therefore M) is 
multiplied by A -1 , the corresponding conjugate field Hq should be multiplied by A in order to keep the 
"scaled" constraint. Since the differential J-^dA of remains locally invariant upon a scale transformation, 
it will be invariant upon all transformations which locally reduce to a scale transformation (plus translations 



and rotations) 119 . This form the class of conformal transformations, which in three dimensions are 
obtained as the group formed by translations, rotations and a three-parameter family of special conformal 
transformations, obtained as the combination of an inversion, a translation, and another inversion. Of course, 
the constraints mentioned above are not invariant in general upon such a transformation. From now on, we 
shall only refer to bilayers, and correspondingly drop the (b) suffix on the rigidity moduli. 
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4 Vesicle shapes 

The theory of the equilibrium shapes of phospholipid vesicles started in 1970, when Canham jl4| showed 
that the characteristic discoidal shape of red blood cells (erythrocytes) can be obtained by the minimization 
of the curvature elastic energy for a particular value of the area and enclosed volume costraints. Actually 
the case of the red blood cells, although of great interest, is somehow misleading, since the spectrin-ankyrin 
network present on the interior of the red blood cell membrane makes its properties somehow different from 
those of a fluid membrane |38], || . 

An extensive study of vesicle shapes was performed by Deuling and Helfrich in the seventies |54[ |33| . The 
problem was reconsidered more recently, in particular because of the experiments of the Sackmann group 
|Q, [7(], [Fill ) which exhibited shape transitions of single vesicles induced by changes in temperature. This 
prompted the research on the identification of the phase diagram of vesicle shapes (for recent reviews, see, 
e.g., |^0[ , and the Proceedings contained in |B2[). The interest has been further enhanced by the observation 
of toroidal vesicles by Mutz and Bensimon 1 92l , which has stimulated several investigations on vescles with 



higher genus [ 107 , 108 



In the spirit of mean field theory, the shape of a vesicle is obtained by minimization of the elastic free 

F = I dA WkH 2 + kk\ . (4.1) 



energy (3.16) 



s 



2 



For phospholipids, k ~ 10~ 19 J |37j, |9lJ| . This expression of the energy also yields a nice argument to 
explain why isolated membranes form vesicles at all. Indeed the elastic free energy F of any given shape is 
independent of its scale: e.g., for a sphere one has 

F = 8itk + 4ttk. (4.2) 

If the membrane is not closed, but has a free edge of length £, the free energy of the edge will be proportional 
to £, times some line tension r, which is of the order of 10 -20 -10~ 19 J nm -1 |33], fl||. As soon as the lateral 
size of the membrane becomes larger than k/t, i.e., 1-10 nm, it becomes energetically favorable for the 
membrane to get rid of its open edge and form a closed vesicle. A detailed study of this transition, analytic 
at zero temperature, and via simulations at T > 0, is contained in ref. |Q. 

The bilayer is relatively permeable to water much less to ions: the permeability ratio is of order 
10 9 . On a time scale of several hours the ions enclosed in the vesicle at its formation remain there. As a 
consequence, any variation in the enclosed volume — due to the permeation of some amount of water — would 
lead to the apperance of an osmotic pressure. Equilibrium is reached when the osmotic pressure balances the 
membrane tension. On the other hand, the time needed to exchange amphiphilic molecules between the two 
layers, or between the bilayer and the solution, are of the order of several hours. Therefore, as long as we 
consider shorter times, we can argue that the enclosed volume V, the total bilayer area A, and the number 
of molecules contained in each bilayer are constant. The last constraint, as we have seen, is equivalent to 
a constraint on the total mean curvature M — § dA H. These constraints can be taken into account by 
Lagrange multipliers, which correspond to the pressure difference p between the interior and the exterior, 
the lateral tension 7, and a quantity /1 proportional to a variation of the spontaneous curvature. In principle, 
one cannot rule out the possibility that a small spontaneous curvature Hq is present. The Euler-Lagrange 
equations read 

S(F + P V + jA + fiM) = 0. (4.3) 

It is a consequence of the Gauss-Bonnet theorem that the Gaussian curvature term does not play any 
role in this variational problem. According to this theorem, the integral ^ dAK is a topological invariant: 

dAK = 4ttxe, (4.4) 

s 

where Xe is the Euler-Poincare characteristic, equal to one minus the number of "handles" of the surface: 

i.e., it is equal to 1 for the sphere, to for the torus, to -1 for the torus with two holes. . . The last term 



in (3.16) remains therefore constant for all continuous deformations of a given surface. The variational 



equation depends only on the parameters p, 7, n and Hq. The solution of the Euler-Lagrange equations will 
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be characterized by the values V of the volume, A of the area, and M of the total mean curvature. If we 
now perform the scale transformation r — > Xr , we have 

V V = X 3 V; (4.5) 

A -> A' = \ 2 A; (4.6) 
M -> M' = AM. (4.7) 

This shape will be a solution of the same Euler-Lagrange equation, corresponding to the same value of k, 
but where 

(4.8) 
(4.9) 
(4.10) 

With these substitutions, the free energy remains invariant. 

We can use this scale invariance to draw the phase diagram of vesicle shape as a function of dimensionless 
variables. One naturally introduces the radius Rq as the radius of the sphere having the same area A: 

Ro =(4z\ ■ (4-H) 



p - 


i 

-> p 


= A- 3 P ; 


7 - 


- i 


= A- 2 7 ; 


M - 




= A-y 



One can thus define: 

• The reduced volume: 

• The reduced total curvature m: 



Air J 



V _ V 

= 6^-t; (4.12) 



M 

m=— . (4.13) 



R 



o 



The resulting Euler-Lagrange equations cannot be analytically solved in general. If one looks for axisym- 
metric shapes one can transform these equations into a system of first-order ordinary differential equations, 
which can be solved numerically |9(|. In this way one obtains the phase diagrams shown in fig. as a 
function of v and Aa = m/Air. The continuous lines (7 pear and C sto denote lines of continuous transitions 
at which the up/down symmetry of the vesicle shape is broken. L poar , L sto and L dumh denote limit shapes. 
The dumbbell region contains for large w-values prolate ellipsoids, and the discocyte region oblate ellipsoids. 

The figure also contains pointed lines, labelled by numbers, which correspond to trajectories observed 
in actual phospholipid vesicles. The photographs of these vesicles are shown in fig. [l8|, together with the 
calculated shapes. In the interpretation of these experiments, the value of the reduced total curvature m 
is inferred by the observed shape changes instead of being measured. In order to explain the transitions 
observed in a given vesicles, one has to assume that the inner and outer layer expand at different rates 
when the temperature is raised, introducing a phenomenological parameter which depends on the vesicle 
one is looking at. The "budding" transition is particularly remarkable, since it implies the existence of a 
"strangled neck" in which the curvature radii go to zero. However, since the curvatures have opposite signs, 
the contribution of the neck to the curvature energy is finite. 

Experimentally, the budding transition appears according to the following scenarios: 

• The vesicle first breaks its up-down symmetry, assuming a "pear" shape, the goes over discontinuously 
to the budded shape; 

• The vesicle goes over continuously via a dumbbell shape to the budded shape. 

While the second possibility agrees with the calculated phase diagram of the BC model, the first one has some 
difficulties: the BC model predicts a continuous budding transition, on the other hand the SC model does 
not predict the intermediate "pear" shape. A model based on area-difference elasticity has been introduced 
to resolve this difficulty [p5[. It seems to me that the continuous elastic model of the membrane breaks 
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Figure 17: Phase diagram of shapes in the BC model. From K. Berndl, J. Kas, R. Lipowsky, E. Sackmann, 
U. Seifert, in L. Peliti (ed.) Biologically Inspired Physics (New York: Plenum, 1991). 



down near the strangled neck, which may have an additional energy cost due to the strain it imposes on the 
phospholipid tails: the extra cost could turn a continuous transition into a discontinuous one. 

Mutz et al. first observed toroidal vesicles in partially polymerized vesicles p2[ . Vesicles of the same or 
higher genus have then been observed in fluid vesicles made with the same phospholipid. This observation 
has aroused great interest, because it opened the possibility of investigating "experimentally" a classic 



mathematical problem posed by Willmore [ 1 1 9 1 . He had asked in fact which shapes minimized, for each 
topological genus, the Willmore functional W = <f s dAH 2 , which is none other than the Helfrich Hamiltonian 
with zero spontaneous curvature. For genus the solution is shown to be the sphere, and for genus 1 Willmore 
proposed the conjecture that the solution was the Clifford torus, i.e., the circular torus in which the ratio of 
the inner circle radius to the outer one is equal to 1/ v2. 

However, all images of the Clifford torus upon conformal transformation correspond to the same value 
of W. There is a family of these images — similar to tori with a non-central hole — known as Dupin eyelids. 
Therefore, the consideration of conformal invariance paves the way to the study of non axisymmetric shapes. 
Upon conformal transformations, the reduced volume v and total curvature m are not invariant. In particular, 
the reduced volume of the Dupin eyelids is always larger that the volume of the Clifford torus, which is equal 
to 0.71. The minimization of the curvature energy finds therefore its solution (if Willmore's conjecture is 
true) if we consider the SC model (without the constraint on the total curvature) and look for surfaces with 
reduced volume larger that 0.71. Other forms can be obtained by the solution of the differential equations, 
as it had been done for vesicles of genus 0. One can check the stability of the axisymmetric forms with 
respect to an infinitesimal conformal transformation which breaks the symmetry. One can thus derive the 
phase diagram shown in fig. [lj] One finds: 

• A region of nonaxisymmetric tori (T na ); 

• A region of discoidal tori (lb;); 

• A region of stomatoidal tori (1st); 

• A region of spheroidal tori (Ts p ). 

The last three zones have boundaries on the limit curves: SToi, where discoidal tori auto- intersect, and Lni, 
-^St, ^Sp where the diameter of axisymmetric tori vanishes. On the continuous lines which separate Tdi from 
Tst and this from Ts p the shapes cross over continuously from one family to the other. All four families 
have the Clifford torus Cliff in common. Circular tori form the boundary between axisymmetric and non 
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Figure 18: Demonstration (1) of a budding transition; (2) symmetric-asymmetric reentrant transition; (3) 
discocyte-stomatocyte transition. Vesicles of DMPC in water: the reconstructed trajectories in (v,m) space 
correspond to the dotted lines in the previous figure. From K. Berndl, J. Kas, R. Lipowsky, E. Sackmann, 
U. Seifert, in L. Peliti (ed.) Biologically Inspired Physics (New York: Plenum, 1991). 
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Figure 19: Phase diagram for vesicles of genus 1 in the BC model. Gray areas correspond to non axisymmetric 
shapes. From F. Julicher, U. Seifert, R. Lipowsky, J. Phys. II France 3 1681 (1993). 

axisymmetric shapes for m > mens = 47r x 2~in^ = 4tt x 1.05. For smaller values of m, this boundary 
is formed by spheroidal tori. Discoidal tori are equilibrium shapes only within the BC model with zero 
spontaneous curvature. 

Experimentally [M one finds indeed that vesicles having the Clifford torus shape turn into non axisym- 
metric shapes (close to Dupin eyelids) upon cooling (which induces a decrease of the bilayer area, and thus 
an increase in reduced volume) (see fig. |20| ). One can also observe axisymmetric discoidal tori, supporting the 
idea that the vesicles are described by the BC model — although one could not rule out the possibility that 
these shapes are merely metastable (see fig. ^T| ). However, one can also observe nonaxisymmetric discoidal 
tori, which do not appear as stable shapes in the phase diagram of the BC model with vanishing spontaneous 
curvature. These shapes cannot be obtained as conformal transforms of a stable axisymmetric shape. This 
indicates that the stability of axisymmetric shapes against infinitesimal conformal transformations is not 
enough to assess their full stability. One can investigate the stability of nonaxisymmetric shapes by taking 
advantage of powerful program Surface Evolver developed by K. A. Brakke 0, which, starting from a 
given surface shape, lets it evolve reducing the Willmore functional at each step. One starts from a shape 
obtained by image analysis of the experiment and looks for the local minimum of the curvature energy 
near to it — and corresponding to the same value of the reduced volume. A nonaxisymmetric discoidal torus 
is obtained, which corresponds to values of v = 0.52 and m/An = 1.11, well in the region of stability of 
discoidal axisymmetric tori. However, its elastic energy is slightly higher and the shape is therefore at most 
metastable. In any case, it will be necessary to take into account the effects of a nonvanishing spontaneous 
curvature in the BC model. 

Vesicles of higher topology have more recently been observed. The important new fact is that for some 
values of the parameters v and m there is a one-parameter family of conformal transformations which 
conserves the Willmore functional and satisfies at the same time the constraints |69t| . 

By numerically minimizing a discretized version of the Helfrich hamiltonian, Hsu, Kusner and Sullivan 
found that the Lawson surface L shown in fig. |22| corresponds to a minimum, with a value F = Fi — 
1.742 x 8itk |p4| . This surface has a threefold symmetry axis and an additional mirror symmetry plane: 
this symmetry is denoted (in the Schonflies notation) by D-ih- Julicher, Seifert and Lipowsky |69|| developed 
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Figure 20: Breaking of axial symmetry of the Clifford torus upon cooling. The transformation is reversible. 
The parameter e denotes the excentricity of the eyelid. The bar corresponds to 10 /im. From X. Michalct, 
Thesis, Paris VII (1994). 




Figure 21: Axisymmetric discoidal tori. In (a) one clearly observes the meridian sections characteristic of 
the discoidal family. Bar: 10 /im. From X. Michalet, Thesis, Paris VII (1994). 
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an algorithm for minimizing the bending energy F for a triangulated surface, and were able to confirm 
the result of rcf. p4j . By applying special conformal transformations to the Lawson surface, one obtains 
a three-parameter family of surfaces (the Willmore surfaces W) having the same value F2 of the bending 
energy. By projecting the three-dimensional space W onto the two-dimensional space (v, to), where v is the 
reduced volume and m is the reduced total curvature, one obtains a two-dimensional region W of the plane 
in which a one-parameter family of special conformal transformations leaves the bending energy invariant 
by still satisfying the constraints. 

Special conformal transformations (SCT) can be parametrized by a vector a = (a Xl a yi a z ), defined by 
the transformation rule 



f/r 2 + a 



{r/r z + a) 

A SCT acting on a surface with initial values v = V\ and to — mi generates a new shape with v = vi(a) and 
to = TOi(a), where 

Ui(a) = v x \l + A^a a + 0{a 2 )] , (4.15) 
mi(a) = toi \l + A^a a + 0(a 2 )] , (4.16) 

in which the coefficients A^ v ' and A^ can be expressed as 

= Q{R A ~R V ) 1 (4.17) 
= 2(R A -R M ), (4.18) 

in terms of the center of volume R v = § dVR/V, the center of area R A = § dAR/A, and the center of mean 
curvature R M = § dAHR/M. Thus, the conformal mode which conserves both v and to can be identified as 
the SCT with a obeying the differential equation 

^=lWxi>», (4.19) 
ds 

where s parametrizes the path in the space W. 

The boundaries of the region W can be determined by first introducing the button surface B shown 
in ^2| This surface is conformally equivalent to the Lawson surface L and corresponds to v = 0.66, and 
to = 1.084 x 47r. It has three orthogonal symmetry planes, i.e., symmetry T>2h- Choose the (x,y) plane to 
be the midplane of the disk, with the centers of the two holes along the x axis. By applying SCT one can 
define a line Cbl of conformally equivalent surfaces of symmetry C% v which connect B to L. This line is 
defined by a = (0, 0, a z ) with < a z < 3.4/i?o- A further increase in a z breaks the threefold symmetry of 
L, generating the line Cls with C2„-symmetric shapes: for a z = 15.5/i?o, the shape along this approaches a 
sphere with two infinitesimal handles at (v,m) = (l,4-7r). If on the other hand we break the x-z symmetry 
plane of B by a SCT transformation with a = (0,a y ,0), we obtain the line Cbs which also approaches a 
sphere at (v,m) — (l,47r). The three lines Cbl, Cls and Cbs form the boundary of the region W. 

Within this region one should be able to observe fluctuations of the shape of the vesicle among conformal 
transforms: this phenomenon has been named conformal diffusion and has first been observed by Michalet 



and collaborators |86|. Examples of this behavior are shown in fig. |23|. There is of course the delicate 
experimental problem of discriminating between conformal diffusion and ordinary thermal fluctuations. 

Outside of W, the shapes of minimal energy have bending energy F larger than Fi- Stability with respect 
to infinitesimal SCT requires that A^"> be parallel to A^ . This is fulfilled by symmetry if there are at least 
two symmetry planes. Thus one can look for shapes which have at least Ci v symmetry. The calculated phase 
diagram is shown in fig. ^J. For comparison, a blow up of the phase diagram of genus-1 vesicles is shown in 
fig. ^5[ One sees that the diagrams look quite similar: in particular the region W corresponds to the line W 
of Dupin eyelids. 

Vesicles of higher genus have also been observed by Michalet et al. [Q. One can expect a larger number 
of transformations which leave the bending energy invariant. Indeed, a different approach |87j is useful: one 
can consider these deformations as positional fluctuations of necks linking two nearby concentric membranes. 
The strategy is to consider the shape of a neck of radius a linking two square parallel pieces of membrane 
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Figure 22: Willmore surfaces with bending energies F = F2 ~ 1.75 x 8itk. These shapes correspond to 
points at the boundary of the region W: The V^h symmetric Lawson surface L; the T>2h symmetric button 
surface B; and examples of Ci v symmetric shapes along the lines Cbs, Cbl, and Cls, respectively. From F. 
Julichcr, U. Scifert, R. Lipowsky, Phys. Rev. Lett. 71 452 (1993). 
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Figure 23: Example of a vesicle of type BS undergoing conformal diffusion. Bar: 10/im. (a): View from 
"above" , showing the absence of symmetry planes normal to the focal plane. A small rotation leads to view 
(b) and eventually to view (c). The series of three views (d) shows the vesicle at intervals of a few seconds. 
From X. Michalet, Thesis, Paris VII (1994). 
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Figure 24: Phase diagram of genus-2 vesicles. Within the region W, the ground state is conformally degener- 
ate. This region is bounded by the lines Cbs, Cbl, and Cls- The Lawson surface L and the button surface 
B are special point at the boundary of W. Adjacent to W, five different regions exist. The symmetries of the 
shapes within these regions are indicated. From F. Julicher, U. Seifert, R. Lipowsky, Phys. Rev. Lett. 71 
452 (1993). 





Figure 25: Phase diagram for genus-1 vesicles. The line W denotes the Willmore surfaces (Dupin eyelids) 
starting at the Clifford torus CL. Compare with the previous figure. 
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Figure 26: Example of a structure with M = 3 membranes and N = 2 necks connecting them. Periodic 
boundary conditions are understood. From X. Michalet, D. Bensimon, B. Fourcade, Phys. Rev. Lett. 72 168 
(1994). 



of size L x L with periodic boundary conditions. The problem breaks into a inner problem, in which the 
surface can be assimilated to a minimal surface (with vanishing mean curvature), and an outer problem 
which can be solved via an electrostatic analogy. The result is that the necks behave as a gas of free particles 
with a hard core repulsion of range £ c ~ v 'aL. Therefore vesicles with a low density of necks will fluctuate 
freely: only when two necks come nearby will they feel the hard core repulsion. 

The approach can be generalized to the case of M membranes connected by N necks: an example with 
M = 3 and N = 2 is shown in fig. |2^. In the limit where M — > oo and N — > oo one recovers the sponge 
phase. 

We have thus seen that the intimations of the Helfrich hamiltonian on vesicle shapes, even at the mean 
field level, reveal an unsuspected richness. 



5 Shape fluctuations in vesicles 

Red blood cells suspended in solution exhibit a remarkable flicker phenomenon, which is best seen with a 
phase contrast microscope and appears as a shimmering at the junction of the rim and the center of the cells. 
The origin of this phenomenon was debated since its first observation in 1890 |12]] and until the quantitative 
analysis by Brochard and Lennon ]TT[ | proved beyond doubt that it was due to Brownian motion. The 
intensity of the phenomenon is due to the fact that the surface tension vanishes. 

Indeed, let us estimate the amplitude of shape fluctuations in a membrane described by the Helfrich 
hamiltonian. We assume that the equilibrium shape of the membrane is planar, and we take its plane to 
be the (x,y) plane. We represent the shape of the membrane in the Monge form, i.e., by giving the third 
coordinate z, as a function of the other two: 



z = h(x,y). 

We set a — (x, y). The hamiltonian then takes the form: 



(5.1) 



F 



J+-K | V 



V/i 



1 + (Vfc) 



(5.2) 



In this equation V denotes the two-dimensional nabla operator, and the term proportional to the Gaussian 
curvature is understood. We assume that the deformation h from the equilibrium shape is small, along with 
its derivatives: i.e., that both slopes and curvatures are small. We can then expand to lowest order in h and 
in its derivatives, obtaining 



F : 



d 2 a 



(5.3) 
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We assume periodic boundary conditions on a square of side L. The Fourier amplitudes h qi where q = (q x , q y ), 
are defined by 

hq = J ^qe^hia). (5.4) 
In terms of these amplitudes, the Helfrich hamiltonian becomes 

F =^E^ 2 + ^ 4 ] K ■ (5-5) 



2L 2 

9 



By the equipartition theorem we have 



h, 



2 , 

2 



k B T 



= L Z 2 " 4 - (5.6) 
7<^ + nq^ 



We see that, if the surface tension vanishes, the fluctuation amplitude diverges, at small q, like q~ 4 . The 
finite size of the membrane imposes a low q cutoff at ir/L. We can then integrate the former expression to 
get an estimate of the fluctuation amplitude: 



2 \ , 2 k B T 



The square amplitude diverges like L 2 . If the surface tension had not vanished, we would have had a much 
weaker (logarithmic) divergence. As an order of magnitude, since n ~ 10~ 19 -10 -20 J, taking L ~ 8/im as for 
human red blood cells, we obtain \/{h 2 ) ~ 0.4-0.7/im. 

More generally, let us consider the fluctuations of a D-dimensional surface in the Monge representation, 
governed by a generic Hamiltonian of the form 

Koz J ^a^hf . (5.8) 

For the case of the rigidity-dominated membranes we have u = 2, whereas for the case of an ordinary 
interface with a nonvanishing surface tension we would have u> = 1. The wandering exponent £ describes the 
behavior of the excursions of the membrane as a function of its lateral size (/i 2 ) oc L 2| » or, equivalently, the 
behavior of the height-height correlation function: 

>(0) -h(a)) 2 ) oc |a| 2C . (5.9) 

We obtain 

C=\{^-D). (5.10) 

We obtain therefore £ — 1 for the case of membranes (whereas £ = for interfaces with nonvanishing surface 
tension), which implies that the aspect ratio of the fluctuations over lateral size does not decrease with 
increasing size. This means that membranes should appear "wobbled" in the same way even at the largest 
scales. As we shall see, the situation is even worse. 

Helfrich |5^] recognized the important effects of membrane fluctuations. O n the one hand, the projected 
area A v of the membrane is reduced with respect to its true area A (see eq. ( 5.37j ) below). Moreover, the 



wandering of the membrane from its equilibrium configuration implies the existence of a long-range repulsion 
between undulating membranes, as it had already been pointed out in 1978 by Helfrich f36f . To fix one's 
ideas, let us consider a membrane constrained by two parallel, planar walls, set at a distance 2d apart. The 
membrane is on average in the middle, but it collides with either wall from time to time. Collisions are 
separated by a typical length L co \\ ~ d 1 ^ , where £ is the wandering exponent. Therefore the free energy per 
unit area will have a repulsive contribution proportional to the density of the collisions, i.e., to L~ ~ d^ D ^ . 
This "steric repulsion" term decreases like d~ 2 for membranes, i.e., just like the van der Waals attraction 
at short distances. One can thus expect in principles regimes in which van der Waals forces dominate, and 
parallel membranes are bound to each other, or in which steric repulsion dominates, and they repel each 
other. The transition between the two regimes is called unbinding: although it is not described by this 
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Figure 27: Thickness fluctuation profile for various values of the ratio e — [iRq/k. From M. A. Peterson, 
H. Strey, E. Sackmann, J. Phys. II France 2 1273 (1992). 

simple argument, as we shall see later, it has been observed in actual membranes. Finally, he recognized 
that the rigidity modulus k itself is renormalized by the fluctuations, and becomes smaller and smaller for 



A more rigorous analysis of shape fluctuations in vesicles requires one to take into account the fact that the 
equilibrium shape is not planar. The necessary techniques were developed by Peterson [f)7j| . These techniques 
have been recently applied to a reappraisal of the flicker phenomenon in red blood cells by Peterson, Strey, 
and Sackmann |)8|. In general one would like to compute the variation of the elastic free energy F under 
a slight deformation of the vesicle shape, starting from a reference equilibrium shape Sq. It is assumed that 
the deformation satisfies the constraints that define the ensemble. There are many ways to parametrize the 
deformation: the most convenient choice is the normal gauge, in which the deformation is described by the 
distance h of the deformed surface S from the reference one So, measured along the normal to So at each 
point. One then computes the variation of F to second order in the deformation: 



The probability of a given deformation is proportional, in the Gaussian approximation, to exp[— 5 2 F /k^T]. 
The constraints can be taken into account in the following way. The deformation h is splitted into a first- 
order term hi, which satisfies the constraints to first order, and a second-order term /12, which is chosen 
to enforce the constraints to second order. The first-order term h\ parametrizes the deformation, whereas 
the variation of the elastic energy is given by the second-order variation of F times h\ plus the first order 
variation multiplied by the second-order term /12. Equivalently, one may consider adjusting the Lagrange 
multipliers in order to keep the constraints satisfied to second order. 

Once the quadratic variation S 2 F of F is obtained, one diagonalizes it by taking advantage of its symme- 
tries: in the case of the red blood cells, the axial symmetry and the up-down symmetry. It is an extremely 
stringent check on the calculation that all the rigid Euclidean motions emerge as zero energy modes. The 
amplitudes of the fluctuating modes are obtained from the equipartition theorem, and are then summed up 
to obtain the fluctuation profile. Since one cannot rule out in principle the existence of a shear elasticity 
modulus fi in red blood cells, due to the loose ankyrin-spectrin network present on the interior of the mem- 
brane, the formalism requires a slight generalization |)7], [9f|. The fluctuation profile then depends on the 
dimensionless parameter e = [iRq/k. In fig. ^?]the thickness fluctuation profile is shown as a function of this 
parameter. We see that, as the shear elasticity is reduced, the fluctuation maximum moves toward the rim. 
At small values of the shear modulus, this effect is more pronounced in the BC model than in the SC model. 
Experimental results are shown in fig. |2^. They are compatible with a very small value of the shear modulus 



larger and larger membranes |5Sj| . 




(5.11) 
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Figure 28: Flicker amplitude over the diameter of a cell. (Arrows indicate the cell boundary.) The observed 
distribution agrees well with the theoretical flicker amplitudes obtained from the BC model. From M. A. 
Peterson, H. Strey, E. Sackmann, J. Phys. II France 2 1273 (1992). 



fj,, with the predictions of the BC model, and with a value of k = 1.4 x 10~ 19 J. This value is of the order of 
that of pure lipid bilayers fl37| , but a factor two smaller than that of a mixture of a physiological lipid and 
cholesterol |}7], ^9) . It is possible that the discrepancy is due to the presence of small proteins |77|] , since red 
blood cell membranes contain about 50% proteins by weight, covering about 20% of the area. 

In order to be able to go beyond the Gaussian approximation, we now need to pay a closer look at 
fluctuations. The most convenient framework introduces the effective potential T[fo] and the renormalization 
group [ p5| , |75| , p8| , [52] ]. Let us consider an fluctuating membrane, whose instantaneous configuration is 
parametrically represented by r = r(a), r € R d and a G R 2 . The free energy T of the membrane is given in 
terms of the partition function 

T = -fc B TTog f Vfexp {-^} , (5.12) 



where F[r\ is the Helfrich hamiltonian 



F j (1,1 



t + ^-H 2 + k q K 



(5.13) 



If we want to assign the membrane a given average shape Fo(ct), we can introduce a fictitious external field 
A, coupling to the position r, such that (r(<r)) A = f (a). The effective potential T[ro] is the defined as the 
Legendre transform of T with respect to A: 



r[ro] = T- J dAX-fo(a). 



(5.14) 



(A precise definition needs to take into account reparametrization invariance.) 

Euclidean and reparametrization invariances dictate the form of the effective potential, just as they do 
for the Helfrich hamiltonian ]95|]. Therefore 

r[f ] = J dA [t cS + ^-H 2 + n cS K +...], (5.15) 

where A is the area element of the average membrane shape, and H, K its curvatures. In this formula 
we have neglected terms which depend on higher derivatives of the average shape. This expression defines 
the effective parameters r c g, k c s and K c ff. If the equilibrium shape is flat, all curvature terms vanish at 
equilibrium, so that 

r cq = r cff Ap. (5.16) 
Therefore r c ff is the effective frame tension of the membrane. 
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Now, the effective frame tension r e s, must be equal to the ll q 2 coefficient," 7, of the two-point height 
correlation function of the membrane, defined by 



W + 0(g 4 ) ' 



(5.17) 



It is sufficient to recall that the effective potential Y is also the generating functional of the vertices, 
and in particular its second derivative yields the inverse of the propagator. By going over to the Mongc 
representation of r we obtain 

-1 8 2 Y 



5h- q Sh q 



(5.18) 



h=0 



We obtain therefore 7 = T e g. In a similar way, the effective rigidity K e s can be related to the coefficient of 
(V 2 ft.) in the effective potential. 

If we wish to investigate the behavior of the membrane at larger and larger distances, we can apply the 
renormalization group approach. The effective Wilson hamiltonian H e s(s) is obtained by integrating out all 
fluctuations whose wavenumbcrs q are contained in a shell A/s < \q\ < A, where s > 1. Again by Euclidean 
symmetry, H c s(s) can be expanded in the form 



n eS ( s ) 



dA 



Teff(s) 



Kes(s) 



H 2 +K eS (s)K + ... 



We then rescale the lengths in order to bring the upper cutoff back to A: 

f — ► s r, , q 

yielding the renormalizcd Hamiltonian 



s 1 q, 



H(a) = J dA 



(5.19) 



(5.20) 



(5.21) 



Since the "order parameter" in our case is the length h, there is no need of the additional "wavchmction" 
rescaling, which usually appears in critical phenomena. The relation between the effective parameters and 
the renormalized ones is simply 



k(s) = k cS (s); k(s) = Kcff(s); r(s) = s 2 r eff (s). 



(5.22) 



The renormalized parameters are actually computed by integrating the renormalization flow equations, 
which are obtained by considering s = 1 + e, where e is infinitesimal. By the definitions of the effective 
potential Y and of the effective Wilson hamiltonian H c s(s) it is clear that 



lim Heff(s) = r. 



(5.23) 



Now, if we are considering an ensemble in which the total area A of the membrane fluctuates and a frame 
tension r is applied to it, we have to choose the initial condition (70, k , Ko) of the renormalization trajectory 
in such a way that 

lim Teff(s) = t. (5.24) 

s— >oo 

We see therefore that the essential step of the renormalization group calculation is the evaluation of H e e(s), 
where s = 1 + e. This is quite similar to the calculation of the effective potential, only that one restricts 
the wavenumbers q of the fluctuations to be integrated out to the shell A/s < \q\ < A. This calculation can 
be performed perturbatively. One chooses as a reference configuration the equilibrium one r an d expands 
in the small (normal) deformation field h = r — tq. The effective Wilson hamiltonian Ti c s{s) can be then 
expanded in powers of the temperature: 



H eS (s) = £ W eff(*) [kpTf. 



(5.25) 
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Figure 29: Definition of the normal gauge. Given a reference shape represented as ro(a), the deformation is 
represented by the displacement h = r(a) — rb(cr) along the normal to the reference shape. 

If one stops at one loop level, it is sufficient to consider only the second variation 5 2 of the Helfrich hamiltonian 
with respect to the deformation field h = ?— ?q. The result reads 

H c ff(s) = F[r ] - k B T\ogjvhexp ■ (5-26) 

The notation J reminds the constraint on the fluctuations. 

The measure factor T>h has to be considered with some 
care |28|, ||^, EH], and has been discussed in the seminar by Thomas Powers in this School ]103[ . In 
principle it contains two contributions: one (known as the Faddeev-Popov determinant) which takes into 
account the volume occupied by different reparametrizations of the same surface. It is defined and discussed 



in Appendix B. The other (discussed in particular in refs. |L3|, 103]), which takes into account the fact that 
the actual number of degrees of freedom involved in the fluctuating surface is proportional to the number of 
molecules, i.e., to the actual area element dA. It turns out, however, that one can simplify the problem (to 
the level of not worrying at all about these contributions) if one works in the normal gauge. We consider 
momentarily membranes embedded in d-dimensional space. Let the reference shape be parametrically repre- 
sented by r = r (a), where r £ R d . The instantaneous shape of the membrane is represented by r — r(<r). 
In the normal gauge one imposes the condition 

(r(a) - f Q (a)) ■ dif {a) - 0, i = 1,2. (5.27) 

The displacement h(a) = r(g_) — fa (a) has therefore d — 2 independent components. In the normal gauge 
the Faddeev-Popov determinant Afp turns out to be local and to contribute only a (divergent) shift to the 
effective membrane tension 7. 

The calculation of the rcnormalization group flow to one loop is reported in Appendix C. The result is 

the following m m, m, H3l : 

dK(s) 3 k B T 

= -TZ , , _/_w..A2 . ( 5 - 28 ) 



ds ~ 47r 1 + t(s)/kA 2 ' 
ffijs) _ 5 k B T 
ds ~ 6tt 1 + t(s)/kA 2 



(5.29) 



dr(s) , , A 2 , , /k + t(s)/A 2 . 

- 2r(.s) + — fc B Tlog -^i . (5.30) 



ds w 4tt V k B T 

The calculation is valid to one loop, hence it only holds if the terms of order k B T are small. This implies 
either k > k B T or r/A 2 > k B T. Let us remark that these equations imply that the bending rigidity k 
decreases when one goes to larger and larger scales: the membrane becomes more and more crumpled. On 
the other hand, the Gaussian rigidity k increases: although this has no effect on the behavior of an isolated 
membrane (whose topology is fixed) , it is important as one considers ensembles of fluctuating membranes, 
as we shall do in the next section. 

This perturbative result allows us to define three different regions characterizing the fluctuations of fluid 
membranes: 



1. Tension- dominated region: 



t > kA 2 ; 



r > k B TA 2 . 



(5.31) 
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Figure 30: Different behaviors of fluid membranes with respect to fluctuations in the (t, plane. 1: 
Tension-dominated region. 2: Rigidity-dominated region. 3: Fluctuation-dominated region. Also shown are 
the renormalization trajectories (i)-(iii) discussed below. 



In this region the fluctuations are small and governed by tension r ( "drumhead model" ) . The flow 
equations can be approximated by 

8k dR „ dr , „„. 

3 OS = s Ts = 0; s Ys = 2T - (5 ' 32) 

2. Rigidity dominated region: 

k > k B T; t < kA 2 . (5.33) 

In this region fluctuations are small, but are dominated by the rigidity term. One has the following 
approximate flow equations: 

dn 3k B T 

'si = -—■ (5 - 34) 

The renormalization of the surface term can be attributed to an increase of the ratio total area/projected 
area due to thermal fluctuations. Indeed, to leading order in k^T this ratio is given by 

+ <»») 



3. Thermal fluctuations dominated region: 

K < k B T; t < k B TA 2 . (5.38) 

In this domain perturbation theory breaks down: one expects that steric interactions and topology 
changes become important. 

As one considers the same membrane at increasing length scale £, the renormalization group flows will 
carry one from one region to another. One can find different kinds of trajectories. Starting from a given 
value kq > k B T of the rigidity one can have: 

(i) If the surface tension To is large enough, the trajectory remains in the tension-dominated region at all 
scales, and is described by the simple drumhead model like an ordinary interface, with rigidity-induced 
corrections with the fixed rigidity kq. 



3G 



(ii) For smaller surface tensions, one may cross over from the rigidity-dominated region to the tension- 
dominated region: it exists therefore a crossover length £ c such that for £ > £ c one is in the tension- 
dominated regime with some (nonzero) effective tension r c g and an effective rigidity K ff . The crossover 
length £ c will be given by 

4 = (5.39) 

V T cff 

For £ < £ c one is in the rigidity-dominated regime: tension effects can be neglected, and the effective 
rigidity depends on the scale £: 

^) = «o-^log(£), (5.40) 

where £q is the molecular size. 

(iii) If the tension is too small, there is a length scale £ c in which one crosses over to the fluctuation- 
dominated regime and the model breaks down. This length can be defined by the condition that k(£ c ), 



as defined by eq. (5.4C), is of order k^T. This length is known as the persistence length [jf5| and depends 



exponentially on the rigidity kq: 

At larger scales the membrane will be crumpled, with a correlation length for its normals of the order 

of e 

Frangois David and Emmanuel Guitter pTj , |30[ | have considered the model with bending rigidity in the 
large d limit. The results can be described as a phase diagram in the tension-rigidity plane. For sufficiently 
large surface tension r, the system is in a "flat" phase with a finite ratio of the projected area to the total 
area. As one approaches a critical line t c (ko), this ratio vanishes: on the other side of this line the surface 
is crumpled and the system does not exist — unless one takes into account steric interactions and topology 
changes. Along the critical line the frame tension remains finite. However, even before reaching the critical 
line, the system becomes instable with respect to nonhomogeneous fluctuations. As one reaches the transition 
line, the frame tension remains finite. 



The effective model describing the membrane at long distances is the Liouville action model [ 100 1 , in- 
troduced by Polyakov in the context of strings. It describes a Gaussian surface (defined by the vector field 
r), coupled to a fluctuating intrinsic metric gij. We discuss it briefly in Appendix D. It is usually believed 
that this model also applies to membranes embedded in finite-dimensional spaces. However, in this case, 
self-avoiding effect are probably essential in determining the actual behavior of the membrane. Moreover, 
it is generally believed that these membranes should assume configurations characteristic of branched 
polymers. Since this is a nonperturbative effect, it is rather difficult to study it analytically. Cates jl8| argues 
that it is related to an instability of the Liouville field theory with respect to the formation of "spikes" . 

A way to control this instability in closed vesicles is to introduce osmotic pressure. Leibler, Singh and 
Fisher [fr9fl verified by numerical studies that two-dimensional self-avoiding vesicles change continuously from 
a deflated state with the characteristics of branched polymers to an inflated state as the osmotic pressure 
p is increased. This transition has been further studied by real-space renormalization group techniques Q 
and by conformal group techniques Jl7[ , which have allowed to determine its critical exponents in 2-D. 
There have been fewer studies of the corresponding transition in 3-D. In particular, Gompper and Kroll ]50| ] 
and Baumgartner || have numerically shown that the deflated-inflated transition is first-order in 3-D. A 
finite-size scaling analysis of the deflation-inflation transition has been recently performed by Damman et 
al. |||. 

I close this section by pointing out that the fact that the effective rigidity decreases at larger and 
larger length scales can be interpreted in terms of the Mermin- Wagner theorem |83), which states that a 
continuous symmetry cannot be spontaneously broken in two dimensions by a local Hamiltonian. In our case 
the continuous symmetry is the Euclidean one, which would be broken for our two-dimensional membranes 
if they exhibited a flat phase. 
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6 Interacting fluid membranes 

We now consider an ensemble of fluctuating, self-avoiding membranes, following Huse and Leibler p5| . The 
effects of excluded volume and of the topological term k J dA K are essential. We shall first go across the 
phase diagram, and identify the phases which may be present. We shall then examine each of them more 
closely. 

Let us first set K = 0, thus effectively neglecting the topological term but allowing for changes in topology. 

From a physical point of view, we are considering a solution of amphiphilc in water. We assume for the 
time being that the bending rigidity kq is quite large, and neglect momentarily renormalization effects. The 
"bare" surface tension To appearing in the Helfrich hamiltonian is proportional to the chemical potential of 
the amphiphile. It can be controlled by changing the amphiphilc concentration (f>. For To large and positive, 
i.e., for a small concentration </> of the amphiphilc, the membrane breaks down into small isolated vesicles: as 
a consequence, the volume is separated into two components: one inside and one outside of the droplets. We 
can assign "spins" to each component (e.g., "down" inside and "up" outside) and we recognize a spontaneous 
symmetry breaking between the two components. 

In principle, we could imagine to force the spins to be up on one side of the sample, and down on 
the other side. In this case there should be a membrane running across the sample, costing a free energy 
proportional to its section. In this phase, therefore, there is a nonzero surface tension t for the membrane 
at a macroscopic level: Huse and Leibler p5[ call it the tense droplet phase. 

Increasing the amphiphile concentration (reducing To) the volume enclosed by the droplets increases. The 
membranes come closer to one another and are likely to form "necks" in order to increase entropy. Therefore 
the region of minority spins becomes connected over longer and longer distances. Well before reaching a 
volume fraction equal to one half, the connected regions can percolate. We thus have two infinite connected 
regions of unequal size. We have obtained a bicontinuous phase. The symmetry between inside and outside 
is still broken. Therefore there will still be a nonzero macroscopic surface tension r between a prevalently 
"up" and a prevalently "down" region. The phase can therefore be called tense bicontinuous. 

If t is further reduced, the symmetry should eventually be restored. We thus obtain an isotropic random 
phase. This phase is analogous to the disordered phase of a ferromagnet, whereas the tense phase is analogous 
to the ordered one. The spin-up spin-down symmetry can be explicitly broken by introducing a spontaneous 
curvature Hq, or a "magnetic held" \x coupling to the spins. The transition from the tense bicontinuous 
to the random isotropic phase should be critical at the symmetry point (Hq — \i = 0) and belong to the 
Ising universality class. Both phases are also called "sponge" phases in the recent literature. In particular, 
the tense bicountinuous phase is called the asymmetric sponge, whereas the random isotropic is called the 
symmetric sponge. They are usually identified with the shear birifringent L3 phase observed in some ternary 
(and also binary) mixtures. 

The tense bicountinuous phase is "paradoxical" : the broken symmetry is not observable. In principle, 
we could have two identical-looking samples belonging to different phases, and one could notice this fact 
only by observing the interface which forms when the two are put in contact. Something similar happens 
in antiferromagnets: the order parameter of an antiferromagnet has a free choice, but it is not possible to 
distinguish between samples corresponding to different values of it. In principle, if one puts them in contact, 
one should observe an interface: in practice, we can probe its order only indirectly. 

Let us now go all the way to a high concentration of amphiphilc. Having to accommodate a large amount 
of membranes, which refuse to cross and pay energy to bend, the most likely organization is to pack them into 
stacks. The membranes are called lamellae in this context, and this phase is called the lamellar phase. From 
the symmetry point of view it is a (lyotropic) smectic A phase, and it is characterized by quasi-long range 
order in the positions of the lamellae, and long-range order in their orientations. The lamellae lie on parallel 
planes, on average, in different regions of the sample. On the other hand, the concentration-concentration 
correlation function of the amphiphile decays (like a power law) to an average value instead of keeping the 
rippled behavior characteristic of the lamellar organization at short distances. 

The quasi-long range positional order can be destroyed by a defect-unbinding mechanism not unlike that 
underlying the Kosterlitz-Thouless transition. One thus goes over to a nematic phase, characterized by expo- 
nential decay of the concentration-concentration correlations (and therefore macroscopically homogeneous), 
but with long-range orientational order. Microscopically, the difference between the lamellar (smectic) and 
the nematic phase lies in the presence of unbound defects. There may be actually different kinds of defects. 
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Figure 31: Examples of defects in lamellar phases: (a) edges; (b) seams. These defects are assumed to be 
forbidden in the simplified model of interacting fluid membranes. 



Figure 32: An edge dislocation which does not break "spin" ordering in lamellar phases. 

The bilayer may end at a free edge, or three bilayers meet on one line (seams), as shown in fig. In 
both cases, by going around the dislocation, the number of lamellae increases or decreases by one. On the 
other hand one may have an edge dislocation, where one lamella bends on itself, as in fig. |3^, or a more 
complicated screw dislocation, but with a mismatch of two in lamella number. 

These defects may also interfere with the "spin" order of the solvent. In a perfect lamellar phase, it is 
possible to assign to regions occupied by the solvent a spin, such that regions separated by a membrane have 
opposite spin. This is not possible if the membrane possesses free edges or seams. However, if these defects 
remain at a microscopical scale, this is still possible to define an ideal surface which spans the defect, and 
the macroscopical "spin" order is not disrupted. With the proliferation of defects, one reaches the point in 
which there is no macroscopically acceptable assignment of spins. In this case, spin symmetry is restored. 
On the other hand also the orientational order characteristic of the nematic phase will eventually undergo a 
transition (of the Heisenberg universality class) . These considerations make it plausible that at a sufficiently 
high temperature (or rather, small bending rigidity k) the nematic phase will leave the place to an isotropic 
phase, which we can identify with the random, isotropic phase discussed above. 

The frontiers between the different phases should correspond to renormalization group trajectories. In 
the large rigidity region, we have seen that the renormalization group trajectories are given by 

dn 3k B T 

US 47T 

dr ( Ak B T\ , , 

"Ts = T { 2+ ^t)- ^ 

Therefore, defining {3 = k^T/n, we have At ~ s 2 +°(/3) ; while A/3 ~ logs. We thus expect the trajectories 
to have the form 

8ttk " 



TCXeXP ^fc^J- (6 ' 3) 
The resulting phase diagram in the (re -1 , r) plane at 7t = is drawn in fig. 

Let us now consider the effects of the gaussian rigidity k, neglecting for the time being its renormalization. 
From the Gauss-Bonnet theorem we know that the integral over the surface of the Gaussian curvature K is 
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Figure 33: Phase diagram for fluctuating membranes at k = 0. The percolation transition between the 
dilute droplets and the tense bicontinuous phase does not correspond to a thermodynamic singularity. From 
D. Huse, S. Leibler, J. Phys. France 49 605 (1988). 



Figure 34: An example of a plumber's nightmare structure. From D. Huse, S. Leibler, J. Phys. France 49 
605 (1988). 

related to the number n c of connected surface components and to the number of handles: 



Let us keep a large value of the bending rigidity k. If we let k become large and positive, structures with 
many handles will be favored: these structures can be built up with minimal surfaces, in which the mean 
curvature locally vanishes. Therefore, for both k and k large and positive, we expect the formation of regular 
structures, made of minimal or almost minimal surfaces, with a great number of handles. These structures 
are collectively known as "plumber's nightmares" , since they represent a pipework in which one does not 
know whether one is in the interior or in the exterior. An example of such structures is shown in fig. |34|. If 
we go on the other side, towards large negative values of 7t, the formation of disconnected components will be 
favored. This can take place in the presence of a rather large concentration of amphiphilc. One will observe 
a large concentration of droplets, which, because of their steric repulsion, will probably tend to organize in 
a close-packed droplet crystal. 

Summing up these considerations, we obtain the phase diagram in the (r, 7t)-plane, assuming a large 
value of the bending rigidity k, and neglecting for the time being the possible occurrence of free edges and 
seams, and their effects. The diagram is shown in fig. |35|. Along the k = axis we recover the sequence of 
phases we had described before. 





(6.4) 
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Figure 35: Phase diagram of an ensemble of nonintersecting membranes in the (t, K)-plane, for fixed k 3> ksT. 
From D. Huse, S. Leibler, J. Phys. France 49 605 (1988). 



We now discuss each of the phases in detail, and take into account the controversial consequences of 
defects and fluctuations. 

Droplet phase 

As I have mentioned above, amphiphilic vesicles are usually formed by sonication and are not assumed to 
be in equilibrium in usual conditions. However, there have been recently observations of droplet phases 
apparently at equilibrium in binary |l6| as well as in more complicated systems |i2], . This phase 

has been named L4. Experimentally, it appears in a very narrow region of the phase diagram, or with a 
rather special amphiphile in the case of the binary mixture. 

In the equilibrium droplet phase one should observe a large polydispersity of droplet size. If one neglects 
the fluctuations, the elastic energy of a spherical droplet is given by 

E d = 8TT^K+^j, (6.5) 

independent of its size. Because of this scale invariance, logarithmic corrections to the free energy of the 
droplet (which are usually neglected in the thermodynamic limit) become very important. One obtains an 
expression of the form 

Fd{n) = adk^Tlmi + F , (6.6) 

where Fq is a constant, and ad is a universal coefficient. As a consequence, the size distribution of the 
droplets becomes (taking into account the chemical potential, p, of the amphiphile) 

Pd (n) CX e -^(«)-Am)/fc B T K n ~ ad e -«/n* ; (g J) 

where n* — k^T/\p\ is a characteristic size, dictated by the amphiphile concentration. Of course, this form 
can only hold for n > n c , where n c is the minimum number of amphiphilic molecules necessary to make a 
droplet. 

The form (fTfl) was first proposed by Helfrich |5ll . It has been reexamined by Huse and Leibler |65| , 



Simons and Cates [ 110 1 , and more recently by Morse and Milner J90j] . The contributions proportional to 
Inn arise from the renormalization of the elastic constant which we have mentioned before, but also by 
the reduction in the translational and undulational degrees of freedom, which depend on droplet size. The 
authors do not agree on the value of ad- Part of the disagreement arises from a controversy over the 
coefficients, a K and a«, which appear in the renormalization of the bending constants. In the previous 
section, I have reported the most widely accepted results, namely a K = —3 and = 10/3. In terms of these 
coefficients, we may write down the following expressions of ay: 



Helfrich |59J, Simons and Cates 11C|: ad = a K + a— /2 = —4/3; 



Huse and Leibler |65||: = a K + + 2 = 1/3; 

Morse and Milner = a K + + 5/2 = 7/6; 

The behavior of the number density of droplets containing n surfactant molecules is rather different, 
depending on whether ad is negative or positive. If it is negative, pd(n) vanishes at small n, exhibits a peak 
for n ~ n* , and eventually drops exponentially. On the other hand, if ad > 0, pd(n) appears to diverge at 



41 



small n, only to be cutoff (by nonlinear effects we have not explicitly considered) for n < n c . For n > n c , 
Pd{n) decreases, first like a power, then (for n > n*) exponentially 

The free energy density (per unit volume) of the dilute droplet phase can be obtained by considering it 
an ideal-gas mixture of droplets of different sizes: 

/ = k B TY,Pd(n) [In (pd(n)A 3 ) - 1 + F d (n)/k B T] , (6.8) 

71 

where A is the thermal dc Broglie wavelength for an amphiphile molecule. Substituting the equilibrium 
values of Pd(n) in this expression, we obtain 

/ = -k B T (/9amph/n* + Pd) , (6.9) 

where p amp h = npd(n) is the number density of amphiphilic molecules, and pd is the total number density 
of the droplets. The number density of amphiphilic molecules, p a mph is related to the volume fraction of the 
amphiphile, 4> by the obvious relation 4> — Pamph^amph, where w am ph is the volume of a single amphiphilic 
molecule. The total amphiphile density is dominated for any ad < 2 by contributions from droplets of size 
n ~ yielding 

n*^(p am p h A 3 eWy/( 2 --\ (6 . 10) 

This yields n* ~ <j/> , where 

f 3/10 for a d = -4/3; 
C = < 3/5 for a d = 1/3; (6.11) 
[6/5 for a d = 7/6. 

On the other hand, if ad < 1, both pd and / are dominated by droplets of size n ~ n*, yielding 

/OC „^(l-«d)/(2-a d ) j (g 12 ) 

which decreases as a power law as <f> increases. However, if 1 < ad < 2, pd and / are dominated by the 
contributions of small vesicles of size n ~ n c , yielding a free energy density that is essentially independent 
of 0. 

In conclusion, I think that a careful reexamination of the calculation od the free energy of fluctuating 
vesicles, in particular taking into account the constraint of the mean curvature, should be done in order to 
clear the controversy. 

Sponge phases 

The volume fraction $ of the solvent contained inside the droplets can be estimated by 

$ oc ^V /2 p d (n), (6.13) 

n 

and is always of order (n*) 5 / 2 . It therefore increases with increasing amphiphile concentration, until it be- 
comes larger than a percolation threshold $ p , which, for a three-dimensional space, will be smaller than one 
half. For 1/2 > $ > $ p , there will be two infinite connected regions occupied by the solvent, separated by 
the membrane, and of unequal volume. This is the asymmetric sponge phase, also called tense bicontinu- 
ous, because the macroscopic surface tension between two phases which differ in the choice of the minority 
component does not vanish. On the other hand, as the amphiphile concentration increases, the symme- 
try between the two components should eventually be restored, leading to the symmetric sponge (random 



isotropic) phase. A compact review of sponge phases is contained in ref. [105]. 

Several experimental observation lead to the conclusion that the shear birifringent L3 phase, which often 
occurs in amphiphilic solutions exhibiting a lyotropic smectic A phase, is indeed a sponge phase. 

• Conductivity measurements one considers an oil-rich phase, where the solvent is an insulator 

and the membrane is an inverted bilayer containing charges coming from the surfactant countcrions. 
The membrane acts therefore as a conductor. The measured reduced conductivity is a fraction of that 
obtained for the same total amount of pure water plus charges. This is orders of magnitude larger than 
the conductivity observed for a nearby droplet phase. This indicates that the bilayers are connected 
and that continuous water paths are present through the sample. 
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Neutron scattering [102, 44j: The scattering intensity can be divided in two regimes: the large wavevec- 
tor q regime, where all data fall on a universal curve corresponding to the scattering from a flat bilayer, 
and the small q regime, where correlations between pieces of the bilayer form a broad peak. In a lamel- 
lar phase the quasi-Bragg peak due to the long-range order of the bilayer is located at q = 2n/d, where 
d = 5/<f> is the repeat distance of the lamellae (6 is the membrane thickness). For a random distribution 
of connected bilayers, one expects d = "/5/cj), where 7 is a number larger than 1. The peak position 
for the L3 phase shows indeed a 1/0 behavior, with 7 = 1.4-1.6. 



Freeze fracture [113 



Direct images obtained by freeze fracture of the L 3 phase show a sponge-like 
structure, quite similar to that of bicontinuous microemulsions, with the difference that, instead of two 
solvents, there is only one solvent present. 

Light scattering data also yield detailed information on the organization of the sponge phase, and will 
be discussed later. 

We can make more precise our description of the sponge phase by defining a lattice model |l^, , inspired 
by the Talmon-Prager [117] and de Gennes-Taupin ]45 ] model of microemulsions, or in a less explicit but 
more general way, by a Landau-Ginzburg approach [ 104 , 2(], |24|, ||| . We shall dwell on the second approach. 
It starts from the consideration that two order parameters are necessary to describe the sponge phase: one 
is obviously the amphiphile concentration 0, or rather, its deviation from a reference value <p*: 



P = <f>- 



(6.14) 



The value (jf corresponds to the volume fraction at a special point of high symmetry, in fact a double critical 
endpoint (DCE). The second order parameter describes the breaking of the symmetry between the two 
infinite connected components of the space occupied by the solvent, which we arbitrarily label "inside" (I) 
and "outside" (O). Denoting their volume fraction by square brackets, we define 



V 



[i] - [O] 



(6.15) 



The underlying I/O symmetry of the hamiltonian implies that rj can only enter the Landau expansion in 
even powers. 

We can now write down the following expansion of the thermodynamic potential <f> = / — pp, where / is 
the Helmholtz free energy density: 



-pp- 



1 



.4 



1 



1 



2^ + 4' + 2^ + 4" +2" ■ 



(6.16) 



This expansion starts from the point in which both p and rj vanish, which (as we shall see) is a DCE. We 
have used some of the available freedom to rescale the fields and the thermodynamic potential in order to 
get rid of irrelevant coefficients: in this way we are left with p, a, and A. This expression is very similar to 
that obtained in the Landau expansion of the Blume-Emery-Griffiths (BEG) model, which was originally 
invented to describe the superfluid/normal fluid transition in He 3 -He 4 mixtures 0. It is useful to switch (via 
a Legendre transformation) to the (a, A, p) space, which is closer to the experimental situation (in which the 
amphiphile concentration is fixed) and also yields a slightly simpler calculation. In this space the transition 
lines are identified by the usual double tangent construction. We therefore consider the following expression 
of the free energy: 

f(p, n) = \p" + -/ + \tf + ir, 4 + ±prf. (6.17) 

If we fix p, and minimize this expression with respect to 77, we obtain a value r\ — ?7m(p) and a corresponding 
value f(p) of the free energy. The double tangent construction can now be applied to f(p). This construction 
corresponds to the requirement of equality both of the chemical potential and of the osmotic pressure of the 
amphiphile: 



f(pi) 



2i 

dp 

Pi — 

dp 



2i 

dp 

ft \ d f 

J(P2) ~ Pi g~ 



(6.18) 
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Figure 36: Schematic phase diagrams obtained from the minimization of the Landau free energy in the (a, p) 
plane, corresponding to the five types discussed in the text. The transition from symmetric to asymmetric 
sponge can be either continuous (dashed line) or first order (solid line). T is a tricritical point, Tr a triple 
point, CEP a critical endpoint, DCE a double critical endpoint, and C a liquid-gas type critical point. From 
D. Roux, C. Coulon, M. E. Cates, J. Phys. Chem. 96 4174 (1992). 



We obtain: 



where p* = —A Therefore 



(a/2)p 2 +p 4 /4-(p*-p)74, \i p<p *. 

This corresponds to a usual symmetry-breaking transition in 77, located at p = p* . For p > p* , one has 
Vm{p) — 0, corresponding to a symmetric state, whereas for p < p* , one has T) m (p) 7^ 0. The rest of the 



calculation involves solving eq. (3.18) where spinodal lines exist. A good qualitative feeling for the phase 
diagram can be obtained just in plotting the function f(p) and looking for double tangents. 

In total, five types of phase diagram are obtained in the (a, p) plane, depending only on the location 
of p* with respect to the spinodal curves. The phase diagrams are shown in fig. one goes from one to 
another by varying p* . The same phase diagrams are represented in the (a, p) plane in fig. The complete 
structure of the phase diagram can be represented in the (a,p*,p) space as in fig. We now describe in 
turn the five types of phase diagram. 

(a) When p* > 0, there is a critical line (CL) for large values of a (a > a* = 1/2 — p* 2 ), separating the 
symmetric phase (rj m (p) — 0) from the asymmetric one. This line ends in a critical endpoint (CEP) 
where the CL meets a first-order line. This first-order line end at a critical point (CP) at p = 0, 
a = 1/2, which corresponds to a demixtion critical point between two asymmetric sponge phases, with 
different amphiphile concentration. For a < a\, the first order line separates a symmetric from an 
asymmetric phase. 

(b) When p* = 0, the line of continuous transition separating the symmetric from the asymmetric phase 
ends at a double critical endpoint (DCE) at p = 0, a = a c = 1/2. There the shape of the coexistence 
curve goes as a c — a cx (p — p*) 2 on one side to oc (p — p*) 1 ^ 2 on the other side. There is no longer a 
separate ordinary critical point. 

(c) When — 1/4 < p* < 0, the symmetric and asymmetric phase are separated by a second order line for 
a > a* 2 = 1/2 — 3/3* 2 , and by a first order line for a < a?,. These lines meet at a tricritical point. 
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Figure 37: Phase diagrams (as above) in the (a, id) plane. From D. Roux, C. Coulon, M. E. Cates, 
J. Phys. Chem. 96 4174 (1992). 



(d) For —1/2 < p* < —1/4, one also finds a first order line between two symmetric phases. This line ends 
at a critical point at a — 0, p = 0. 

(e) When p* < —1/2, the tricritical point T disappears and the continuous transition between symmetric 
and asymmetric phase ends at a critical endpoint (CEP). 

Upon looking at this phase diagram it is possible to identify the starting point of the expansion as the 
double critical endpoint (DCE) appearing in the diagram of Type (b). This high order critical point has 
in fact an "upper critical dimension" of 8/3. Therefore, the mean field treatment describes correctly the 
behavior around this point for three dimensional systems. On the other hand, the corrections to mean-field 
theory are very weak (logarithmic) around the tricritical point. One should thus expect the mean field 
predictions to be satisfied in a wide region for systems which happen to be near the DCE or the tricritical 
point T. 

The occurrence of a line of tricritical points in the (a, A, p) plane implies that, even in a binary surfactant- 
solvent system there is a finite probability of finding a tricritical point in a (0, T) phase diagram, without 
any special adjustment of other thermodynamic variables. 

Since the Landau expansion obtained above has two independent parameters (up to rescalings), and 
since the elastic description of bilayers has two independent elastic coefficients k and k, it should be possible 
to span the whole of the phase diagram by varying k and k. In practice one can vary these quantities by 
adjusting concentrations of cosurfactants, salt, etc. These systems may therefore provide a testing ground 
for theories of higher order critical behavior. 

In order to investigate the behavior of correlation function, one adds suitable gradient terms to the 



expansion (6.16). They involve the gradients of the two order parameters p and 77, although only the p-p 



correlations are directly observable. The most general form of the gradient terms is 

A$ = |i|Vp| 2 + ^|Vr,| 2 + ^Vr? • Vp|. (6.21) 
In the simple case in w hich 7 p = 7 C = 0, it i s pos sible to calculate exactly the correlations within the 



Gaussian approach [104]. The third term in eq. (6.21) is not quadratic in the fluctuating fields and prevents 



the use of the Gaussian approach when j c ^= 0. In this case, the probability distribution for p(r) depends 
only on the local value of T](f). One can then integrate over the p(r) field, obtaining an effective potential 
for the r\ field: 

r[u]= /dfj^ + ^ + lMVr,! 2 ], (6-22) 
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Figure 38: Representation of the Landau phase diagram in the (a, /x, A) space. Continuous line: line of 
critical points. Broken line: line of tricritical points. Thick line: line of triple points. Dotted line: line of 
critical endpoints. DCE: double critical endpoint; CE: critical endpoint; TCE: tricritical endpoint. From D. 
Roux, C. Coulon, M. E. Cates, J. Phys. Chem. 96 4174 (1992). 



where A' = A + fi/a, and A = 1 — a/2. This form implies 



(7/(^(0)) 



4w 'j n r ' 



where £ r) = Jrj/A'. 

The amphiphile concentration-concentration correlation function is defined by 

g(r) = {5p(f)5p(0)) , 



(6.23) 



(6.24) 



where 5p(r) = p(r) — (p) . It contains two contributions: the direct contribution, which in our hypotheses is 
very short-ranged (delta-like), and an indirect contribution, due to its coupling to the r\ fluctuations. It is 
not difficult to see that one has in fact 



9(f) = W) + JL [(^(rVCO)) 2 - (V 2 )' 
Since we assume gaussian fluctuations, we have 

(V 2 (r)v 2 (0)) = (v 2 ) 2 + 2(v(f)v(0)) 2 - 
By using the expression of the correlations of 77 obtained above, we have 



9(r) = 

a 



) 1 e - 2r /«" 



whose Fourier transform is 



where 



I{q) = Ci 



C 2 



4a 2 (4irj v r) 2 ' 
tan- 1 (^, ) /2) 



C 



C 2 



167T07 2 



(6.25) 
(6.26) 

(6.27) 

(6.28) 
(6.29) 



This expression can be compared with the results of light scattering experiments. It involves a 1 /g+const. 
behavior at large q instead of the more usual (Ornstein-Zernike) one: l/q 2 + const. The result should hold 
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when the hypothesis of neglecting 7 P and 7 C applies: i.e., near the critical symmetric-asymmetric line and 
away from the tricritical point. 



A perturbation method 1 105, App. A] can be introduced to obtain the expression for I(q) in the general 
case. The result reads 

1{q) = {M = + K ( a+lpq >r ^7^' (6 ' 30) 

for the symmetric sponge phase, and 

1 



(l + ( 7c9 2 /2))' 



(a + 7p g 2 ) 



2 \2 



KAv) 2 , „ tan- 1 ( 9 ^/2) 

A2- 



(6.31) 



for the asymmetric sponge phase. 

Two different experimental situations have been analyzed by light scattering. In the first case |104|, 
one has studied an oil-rich sponge phase as a function of membrane concentration (this experiment will 
be denoted as "oil" dilution). In the system under investigation, the sponge phase undergoes a first-order 
transition to a phase of practically pure solvent. In the second case [^4], [25|, one has studied a water-rich 
system which exhibits a second-order phase transition between a symmetric and an asymmetric sponge 
phase. The transition is made apparent by a maximum turbidity line (MTL) separating the two isotropic 
phases. This is the line where the structure factor I(q) (proportional to the osmotic compressibility) diverges 
at q = 0, corresponding to very large fluctuations in the amphiphilc concentration. In this way the rather 
remarkable feat of exhibiting a symmetry breaking which cannot be directly observed has been achieved. 
This experiment will be denoted as "water" dilution. More recently, there have been careful experiments on 
the quaternary system Ma-octylbenzenesulfonate (OBS)/n-pentanol/brine fjofl . These experiments challenge 
the picture of the asymmetric-symmetric transition given above, and will be discussed later. 

The inverse scattering intensity 1/I(q) measured in oil dilution experiments is shown in fig. || as a 
function of the wavenumber q. The form is strikingly different from the Ornstein-Zernike (dotted line), and 



fits well to the simple form ( 6.28 ) (solid lines). The fact that the simple expression works so well indicates 



that the correlation length for p is rather small. It is possible to extract from the data the ratio of the arctan 



to the constant term appearing in eq. (6.28). This ratio turns out to be about 2/3, and is a measure of the 
relative importance of the fluctuations of rj with respect to the fluctuations of amphiphilc concentration. 

The osmotic compressibility 1(0) as a function of the amphiphilc volume fraction (f> varies approximately 
as l/</>, but with logarithmic corrections: 

7(O)c<(01og(^*)). (6.32) 



This behavior can be understood in terms of the scaling laws ]6f| 101 of the free energy, which derive 
from the scale invariance of the Hclfrich hamiltonian and the renormalization of the elastic constants. As 
we have seen, if all lengths are dilated by a factor A, the Hclfrich hamiltonian remains invariant. In this 
transformation, the volume ratio <fi remains invariant, although the surface area per unit volume is decreased 
by a factor A -1 . However, the membrane thickness <5 is also scaled by A. Microstates of the original and of 
the rescaled system have the same energy, and therefore the two systems should possess the same free energy. 
In order to use this argument to describe the behavior of the same system at different dilutions, we observe 
that the dilated system (with volume ratio (j)) corresponds to the diluted system (with volume ratio </>/A) 
except that its thickness S has been rescaled. Now the thickness appears as a cutoff in the renormalization 
of the elastic constants, which has we have seen produce logarithmic terms in the free energy. As a result, 
we expect the free energy F of the solution to depend on the volume ratio </> in the following form: 

F = k B T(/) 3 (A + B log (/)), (6.33) 

where the constants A and B depend on dimensionlcss combinations of the elastic constants and of k^T. 



This implies the behavior (3.32) of the osmotic compressibility. The logarithmic corrections observed in 1(0) 



are therefore an experimental evidence for the renormalization of the clastic constants. 
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Figure 39: Plot of the inverse scattering intensity against q for a sponge phase at different concentrations. 
Dotted line: Ornstein-Zernike form. Solid lines: cq. (133). From D. Roux, C. Coulon, M. E. Cates, 
J. Phys. Chem. 96 4174 (1992). 



The water dilution experiments allow to exhibit the divergence of the correlation length £ v and of the 
osmotic compressibility 1(0) at the symmetric-asymmetric transition. Experimentally in particular 1(0) oc 
(/i — /ic) - *", where C — 0-5. The Landau-Ginzburg approach predicts 1(0) oc (fj, — pc)^ 1 ^ 2 . However, the 
relation between density and chemical potential is singular at the transition, and the exponents undergo 
the so-called Fisher renormalization (4lj. One has 1(0) oc (p — Pc)^ 1 near the tricritical point. Near the 
critical point one should use nonclassical (Ising) exponents, and the Fisher renormalization has only a small 
effect. The result reads 1(0) oc (p — p,,) - "/^ 1 ""', where a ~ 0.11 is the specific heat exponent. This yields 
£ ~ 0.12, also in contrast with the measured value. This discrepancy may arise from a crossover from critical 
to tricritcal behavior, but has not yet been solved. 

Lamellar phases 

Membranes can assume two forms that minimize the bending energy: as a stack of planar sheets, or as 
periodic structures of minimal surfaces (plumber's nightmare). In both cases the mean curvature vanishes, 
and therefore the free energy of the structure is determined by the contribution of the gaussian rigidity and 
by fluctuation effects. Lamellar phases have been observed for a long time, and they are usually known as 
L a phases p4j . 

The stability of the lamellar phases arises from a compromise between the area coefficient To, which wants 
to pack the layers as close as possible, and the interaction between the lamellae. The interaction between 
neutral membranes (for the electrostatic properties of membranes I refer to the review by D. Andclman jlj]) 
arises from the following sources: 

(a) A strong, short range repulsion due to the hydration interaction, yielding the interaction potential per 
unit area (as a function of the membrane spacing d) V\ L (d) — exp(— d/X^). 

(b) A long range attraction due to the van der Waals forces: 

2 1 
(d + S) 2 + (d + 2S) 2 \ ' 

for d > 5. (6.34) 
Here 5 is the bilayer thickness and the parameter W is known as the Hamaker constant. 
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(c) The long range steric repulsion, due to the occasional collision between lamellae f5f|: 

v *w = *w=$- (6 - 35) 

Here ch is a universal coefficient, over whose value there is some controversy J5(| |29|| . 

However, one cannot simply sum up these contributions to obtain an effective potential V c g (d). This potential 
has a minimum either at very short distances or at infinite separation: as a consequence one would have a first 
order transition between lamellar phases at high amphiphile concentration and a practically pure solvent. 
On the other hand very dilute lamellar phases (where d can be as large as 1 m/i) have been observed j76|. 

Lipowsky and Leibler |8l| introducend a complicated functional renormalization scheme to attack the 
problem, and obtained a continuous transition between a state of "bound" lamellae and one of "unbound" 
lamellae, as the Hamaker constant W crossed a critical value W c . In this theory, the mean lamellar spacing 
d diverges like a power law near the transition: d ~ (W — W c )~^ , where ip = 1.00 ± 0.03. This result did not 
appear to admit a simpler description. 

Recently, Milner and Roux |38| as well as Helfrich [[30| have proposed simpler (Flory-like) explanations 
of this continuous "unbinding" transition. However, the theory of Milner and Roux applies to the uniform 
lamellar phase, whereas the one by Helfrich describes the unbinding of a finite stack of n layers from a 
wall. We shall thus dwell on the theory by Milner and Roux. First of all, one chooses the volume ratio cf> 
of the amphiphile as the order parameter. Then the free energy is obtained in analogy with the van der 
Waals theory of the liquid-gas transition, by adding to the Helfrich estimate of the entropy of the stack 
(which only takes into account the hard-wall repulsion) a term representing the correction to the second 
virial coefficient due to the other interactions. In analogy with the Flory-Huggins theory of polymers, the 
entropy term is estimated in the absence of all correlations due to the attractive interactions, and in the 
correction to the virial coefficient one neglects the connectivity of the surface. This enthalpic correction is 
referred to "patches" of the membrane of size v = Iq8, where £q = (n/k-oT)8 2 is the shortest length where 
the membrane could bend. The result reads 

/lamW = '"^V ^ k B T X ^. (6.36) 

The coefficient x appearing in the enthalpic correction can be given as a function of the interaction U v {f) 
between patches by 

1 

2V 2 

as long as the interaction between patches is sufficiently weak. We expect that \ vanishes linearly as the 
Hamaker constant reaches some critical value. The phase diagram can be obtained by minimizing with 
respect to <j> the free energy 

g(x, /«) = /iam(</») ~ (6-38) 
where /i is the chemical potential of the amphiphile. The resulting phase diagram is shown in fig. The 



X = -— I dV {1 - cxp [-U„(r)/k B T]} , (6.37) 



expression (3.3£) has the form of a Landau expansion around <f) — 0. In this expression, the cubic term in 
4> does not necessarily lead to a first order phase transition, because negative values of <f> are unphysical (a 
similar situation arises in the mean field treatment of percolation) . Therefore the order parameter vanishes 
at the second order transition with an exponent 1 instead of the usual mean field exponent 1/2. 
Topological instabilities 

Wc have thus found that it is possible to stabilize lamellar phases with arbitrarily small amphiphile con- 
centrations, at least if the Hamaker constant is large enough. However, we have argued before that as the 
solution is diluted one should go over to isotropic or cubic phases. In order to understand this transition 
it is necessary to include topology changes and renormalization effects, which are absent from the previous 
considerations. Two recent papers deal with this problem [E8| |8^]. The essential point is that, while the 
bending rigidity k decreases as one goes to longer and longer length scales, the gaussian rigidity H increases. 
Therefore structures of higher and higher genus are favoured at longer and longer scales. 

In ref. f89fl one starts from the expression of the Helfrich hamiltonian in terms of the principal curvatures 
Ci and C2'. 



F = dA 



(Cl + C 2 ) 2 + KClC 2 



(6.39) 
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Figure 40: Phase diagram of a lamellar phase in the (/x, x) plane. Broken line: first order transition. Solid 
line: second order transition. Thin lines: spinodals. From S. T. Milner, D. Roux, J. Phys. I France 2 1741 
(1992). 



and rewrites it introducing the topological rigidities n + = n + k/2, k_ = — k/2: 



F = dA 



(d + c 2 f + -K- (ci - c 2 ) 2 



(6.40) 



One has of course k = n + + When either topological rigidity vanishes, the simply connected membrane 
becomes unstable with respect to fluctuations leading to a different topology: 
For k+ < 0, towards the formation of many spherical surfaces; 
for K- < 0, towards the formation of an infinite minimal surface. 

Since k = n + + k_ , this instability occurs before the instability towards a continuous deformation of the 
membrane, which occurs when k becomes negative, except in the special case in which n + ~ k_ ~ 0. 
The topological rigidities renormalize with respect to the scale I according to the equation 

« ± W = K± -^lo g y . (6.41) 

In this equation, £ is a microscopic distance (of the order of the lateral size of an amphiphilic molecule), k± 
are the bare values of the topological rigidites, and the coefficients a + = 4/3, a_ = 5/3 can be expressed 
in terms of the renormalization coefficients a K = —3, a— = 10/3 for k and k. This renormalization group 
equation defines two characteristic lengths £± oc exp (Attk±/ (a±k-QT)) for the onset of the topological in- 
stabilities. One expects that a lamellar phase will melt when the lamellar distance (or, better, the repeat 
distance d = 8/<j>) becomes of the order of the smaller of these lengths. 

In re f. |4q | a more detailed physical picture of this transition is contained. One can distinguish three 
regimes 



Regime A 



k < -q-K, i.e., £ K <C £ K ; 

10, 



Regime B: -7t > fK, but |k + fn\ > 0(1) k B T, i.e., £ K < 
Regime C: \k + f \k < 0(1) k B T, i.e., ^ ~ 6f- 

The author considers the elastic energy of topological defects. These defects can be in particular passages 
connecting two lamellae, or droplets confined between the lamellae. 

Let us first consider a passage connecting two lamellae a distance d apart. Its curvature energy can be 
estimated as 

E p =0(1)k(£)Q 2 ~47r«(r), (6.42) 
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where r is the size of the passage neck, and £ is the lateral size of the passage deformation region. At distances 
from the neck smaller than £, the presence of the passage curves the membrane giving rise to the first term in 
eq. ( 6.42| ). The second term in this equation is the gaussian curvature contribution, and originates essentially 



from the region of the neck. Therefore the two rigidities are evaluated at different scales. (This effect is not 
taken into account in rcf. [B9J ) . The passage deformation region is approximately a catcnoid, yielding the 



condition d ~ 2r\og(£/r), relating r and I. We let this constraint in eq. (6.42) and obtain the equilibrium 
energy and size of the passage, as a function of d. Assuming n(d) 3> k B T, and neglecting slowly varying 
logarithmic factors, we obtain 

E p = -a-fc B Tlog^-, (6.43) 

4 q = 0(1) (gf) 1 ^ d, (6.44) 
r eq = 0{l)d. (6.45) 

The elastic energy of the passage decreases as d increases, and becomes negative as d becomes of order 
On the other hand, the deformation region £ eq increases with d. 

The interaction among passages can be represented by a hard core repulsion, with effective hard core 
size £ cq . When the area density of n p of passages is smaller than £~^, a simple grand canonical calculation 
yields n p {d) ~ no(d), with 

which increases with increasing d. Such a growth of the passage density was observed by Harbich et al. in 
egg-lecithin membrane stacks [p4j . The free energy of the "passage fluid" per unit membrane area is given 
by F p (d) = —kBTn (d) ~ — d 4 ^7The passages begin to overlap when d becomes of order d* , where 

7 T \ 3/ Q 7T 

knT 



d '=m>) (6 - 47) 

In this situation, the free energy can be estimated by minimizing with respect to the passage density n a 
trial free energy containing the contribution of the hard core repulsion: 

F p (n,d) = k B Tnlog ( —^j^] - k B Tn log (l - n£ 2 c ) . (6.48) 



The first term reproduces the dilute result n — no(d), while the second represents the free energy increase 

~ £~ 2 



due to the hard-core repulsion. We obtain therefore n p ~ £ c 2 for c? 3> d* , with a free energy 



F p (d) ~ -O(l)(^log(no(d)^ 



"cq; 



0(l)M^log(l + (d/d*)^)- (6-49) 



This free energy corresponds to an attractive interaction between membrane pairs, which tends to keep them 
at a distance d ~ f^. This attraction actually dominates over the Helfrich steric repulsion 

V st (d) = -ch^J^. (6.50) 

for d > 

The free energy of the lamellar phase per unit volume is given by 



-Flam (to) = min$(d,ro), 

d 



(6.51) 
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where tq is the area coefficient (proportional to the chemical potential of surfactant) and 

$(d, T ) = i [r + F p (d) + V st (d)] . (6.52) 
a 

This expression does not include the van der Waals interaction. In the absence of the passage term, one has 
a single minimum for To < 0, at d — d cq ~ |to| -1 / 2 . One would thus expect a swelling of the lamellar phase 
until d ~ £ K , where the lamellae crumple. But the inclusion of the passage interaction changes the behavior, 
and d cq remains essentially fixed at also for positive t , where one would expect a diluted droplet phase. 
The free energy of the droplet phase can be shown to be essentially zero compared to that of eq. ( J6.51 ) . The 



lamellar phase therefore gives way to the droplet phase when the free energy ( |6 .51 ) vanishes, what happens 



for t ~ (fceT) 2 / 6c) 2 • This only holds in regime A, where £- <C £ K . In regime A one therefore expects 
the lamellar phase (with a repeat distance d ~ d*) to remain stable also for a positive area coefficient. One 
would thus have a strong first-order phase transition between a passage-rich, but still oricntationally ordered, 
lamellar phase and a dilute droplet phase. 

In regime B, on the other hand, passages will be very dilute and the most important role is played by the 
droplets. The length scale of droplets is fixed by , which is much smaller than £ K in this regime. In ref. |4f| it 
is assumed that the elastic energy of a droplet of scale R is given by E d (R) = 8ttk + (R) = 2a + k^T \og(^ + / R) . 
(The results would probably have to be modified in view of the results of ref. |MJ.) By means of this 
expression, it is possible to compute the free energy, per unit area of lamella, of a polydisperse solution of 
droplets confined between two lamellae. The result reads 

F d (d) = dF d + f d (d), (6.53) 

where T d is the free energy per unit volume of the droplet phase, which in our hypotheses is given by 



The effects of the confinement are contained in f d (d), which is given by 

U(d) = 0(1)^ (£) 2Q+ A' 2 ~ d 2/3 , (6.55) 

corresponding to an attractive interaction between the lamellae. Therefore the difference between the free 
energy density of the lamellar and the droplet phase exhibits a minimum at a value d* = d*(jo). For a critical 
value t c ~ — (fcBT) 2 /(K(d max )(i 2 nax ) of To the free energy difference vanishes and the lamellar phase melts 
into the droplet phase. The maximal distance d max between the lamellae is given by 



3/ a 



+ 



<U« = d*(T c )~ ?d~£d- (6.56) 

This is also the radius i? max of the largest droplets. 

In this regime it should therefore be possible to observe a weak first-order transition to a dense droplet 
phase. This transition takes place when d reaches the value £ + . A transition to an isotropic sponge phase 
can only occur in regime C (corresponding to d max ~ -R m ax ~ ~ ~ £+) where the creation of passages 
between the largest droplets is favored. 



These expectations agree with those of ref. [89 in regime B and regime C. The predictions for regime A 
are different. In particular they contradict the prediction made in ref. p9| of a transition from the lamellar 
phase into a L3 phase of multiply-connected surfaces with small mean curvature. 

Free edges and gauge theories 

I discuss in this subsection the elegant connection made by Huse and 



Leibler |66 between the sponge phases of membranes and the Z2 gauge- Higgs models on a lattice. 

When one considers real membranes, one should also take into account the possibility of free edges. Free 
edges (see fig. ^T](a)) involve a region with strong curvature in the amphiphilic arrangement, and where 
contact between water and the tails of the amphiphilc is cased. It is likely therefore that these defects cost a 
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Figure 41: Theoretical construction for the definition of the macroscopic line tension. An edge (thick line) 
is forced in the sample, connecting points B and C along the curve C\. The return path C2 may either be 
a straight line (a) or a curved line close to C\, depending on the energy cost of the creation of membrane 
patches. From D. H. Huse and S. Leibler, Phys. Rev. Lett. 66 437 (1991). 

relatively high price in terms of energy per unit length. However, this price is not infinite, and it is interesting 
to ask how the phase diagram of the membranes is affected by introducing them. 
One can stard by adding to the Helfrich hamiltonian a term 




(6.57) 



where the integral is over all edges, and Ao is the chemical potential per unit length, of the edges "living" 
on the membrane. For infinite Ao edges are suppressed and the previous results hold. When Ao is finite, 
we start running into difficulties in separating without ambiguities "inside" (I) from "outside" (O) in the 
sponge and droplet phases. Moreover, it is no more possible to define an effective surface tension, as we had 
done in the gedankenexperiment mentioned in the discussion of the asymmetric sponge phase. Indeed, if we 
try to force a macroscopic interface between a prevalently "up" and a prevalently "down" spin phase, the 
system will prefer to introduce a macroscopic hole in the membrane, and the free energy will be higher by 
a term of order L (instead of order L 2 ) with respect to a homogeneous system. 

Instead, one may consider the effects of forcing a macroscopic edge through the system. Imagine to 
introduce an edge along a curved line C\ connecting points B and C, set a distance I apart, as in fig. ^l|. The 
length of the curved line is equal to L. The system must choose the return path C2 connecting C to B. If the 
free energy cost needed to create a macroscopic patch of the membrane is small, nothing prevents the system 
from choosing the straight line BC as C2 (fig- [lT](a)). In this case, the free energy price of the free edge BC is 
proportional to £, and the macroscopic line tension is nonzero. On the other hand, if the free energy cost of 
the membrane is large, the system prefers to return along the curve C2 situated very close to Ci (fig. 0( b ))» 
in order to avoid creating extra surface. In this regime, the macroscopic line tension is formally zero. The 
first regime is expected to hold in the symmetric sponge phase. We expect the macroscopic line tension to 
vanish in the asymmetric sponge phase and also in the phase corresponding to the droplet phase, but in 
which small membrane patches with free edges are also allowed. This is the droplet- and- disc phase. 

In order to understand more closely the phase diagram, let us fix the value of the microscopic line tension 
Ao, of the rigidity Ko, and let us start from a large value of the area coefficient To. We shall be in the droplet- 
and-disc phase, and the edges appear as closed rings, spanned by the membrane. As To decreases, an infinite 
piece of connected surface appears, although (with sufficiently large Ao) the size of the edge loops remain 
finite: this is the asymmetric sponge regime. If we now reduce Ao, one reaches a regime in which connected 
edge paths of infinite length are present: this is the sponge-with-free edges regime. One can picture this 
regime as a melt of ring polymers (the edges) of arbitrary length. These polymers are Gaussian at large 
lengths, and therefore have only a finite probability of coming close to themselves. However, there are an 
infinite number of these polymers, and different polymers will come close to each other infinitely many times. 
So the surface attached to one edge will almost certainly be connected to the surface attached to the other 
edge, and therefore one expects only one connected piece of membrane to be present. 

We have therefore four different isotropic regimes: droplet-and-disc, asymmetric sponge, sponge with 
free edges (all three of them with a vanishing macroscopic line tension), and the symmetric sponge regime 
(with a nonzero line tension). Only the symmetric sponge regime is separated from the others by a bona 
Gde phase transition. 
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Figure 42: Schematic phase diagram for an ensemble of membranes with edges, according to D. H. Huse 
and S. Leibler, Phys. Rev. Lett. 66 437 (1991). The dashed lines correspond to geometric (percolation) 
transitions, while the solid lines show the thermodynamic phase transitions. The tq > of the phase diagram 
corresponds to the phase diagram of the gauge Higgs model shown in the inset. In the inset the dotted 
line is the self-duality line; the unconfined-confined transition on the vertical line has its dual equivalent in 
the unconfined-Higgs transition on the horizontal line. 



For sufficiently large and negative tq one goes over to the orientationally ordered phases. Here, the 
presence of free edges is incompatible with smectic (positional) order, as we have discussed before: but we 
may have, or not have, a macroscopic line tension. We have therefore two nematic phases, depending on the 
value of the macroscopic line tension: nematic with line tension, and nematic with free edges. The resulting 
phase diagram in the (to, Ao) plane is sketched in fig. 

The part of the phase diagram involving the sponge phases is in fact equivalent to the phase diagram of 
the Z 2 gauge-Higgs system (for a review, see, e.g., p5[). The connection can be made explicit by considering 
random surfaces on a lattice. The lattice length should correspond to the persistence length £ K : at longer 
lengths, the bending rigidity can be neglected. Thus, let us consider a random surface living on the plaquettes 
of a simple cubic lattice. Label the elementary cubes by indices i. One introduces Ising spins <7j = ±1 on 
each elementary cube, and a gauge field f/y = ±1 on each plaquette ij (separating cube i and cube j). 

The plaquette ij is occupied by the membrane if diUijUj — — 1. The link shared by cubes i, j, k, and I 
is occupied by an edge (or a seam) if the number of adjacent plaquettes containing membranes is odd. This 
corresponds to UijUjkUkiUu = — 1. The Hamiltonian corresponding to these fields is therefore 

H = -tq OjUjjaj - A ^2 U lJ U ]k UkiU li . (6.58) 

(ij) (ijkl) 

This hamiltonian is invariant under the gauge transformation 

Cj — > e^; Uij — > ei€j Uij; e t = ±1. (6.59) 

Because of this gauge invariance, each configuration of the membrane corresponds to 2 N configurations of 
the lattice gauge model. For infinite Ao, the gauge may be fixed so that all the gauge fields C/y are equal 
to +1, and one recovers the usual Ising model for the spins (7j. This has a critical point as a function of to, 
corresponding to the symmetric-asymmetric sponge transition. 



However, the model (6.59) is self-dual and the Ising line at infinite Ao is mapped to the pure gauge line 
To = by duality. There is correspondingly a phase transition from the symmetric sponge to the sponge 
with free edges, also belonging to the Ising universality class. 
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Away from the tq = or Ao = oo lines one may exploit Monte-Carlo simulations of the Z2 gauge 
system p8[ . The result is that these lines of continuous transition become first-order line near the triple 
point where the three phases meet. A line of first-order transition also extends along the self-dual line, 
separating for a while the asymmetric sponge from the sponge with free edges phase, and terminating in a 
critical point, beyond which there is no thermodynamic singularity separating the two regimes. 

As I have mentioned above, the experimental results of ref. pOfl on the quaternary OBS/?i-pentanol/brine 
system could not be interpreted without diffi culties in the framework of the phenomenological theory of the 
asymmetric/symmetric transition of ref. |l05| . In this system, one observes a sharp maximum of the osmotic 
compressibility 1(0) on the maximum turbidity line (MTL) far inside a one-phase region; however, the light 
scattering function I(q) has a Ornstein-Zernikc q~ 2 form on the n-pentanol rich side of the MTL (while it 
deviates from this form on the n-pentanol poor side). 

This behavior is not compatible with the predictions of ref. |105| , which would imply a (weak) divergence 
of 1(0), and a different form of the scattering function, nor with the hypothesis of a tricritical behavior, 
because one is deep inside a one phase region (while a tricritical point is necessarily close to a two-phase 
region). The authors suggest to introduce a different phenomenological model, introducing a new order 
parameter, to, linearly (and not quadratically) coupled to the amphiphile density p. In this way one would 
explain both the strong divergence of 7(0) and the Ornstein-Zcrnike form of I(q). The divergence is then 
"rounded off" because the amphiphile concentration acts as a "magnetic field" for m, and the critical point 
is reached only for that particular value of the concentration for which the "magnetic field" vanishes. 

The authors then go on to suggest that the MTL in this system can be related to the symmetric 
sponge/sponge- with- free-edges (S/SFE) transition considered by Huse and Leibler Q. The idea is that 
the (small) alcohol molecules essentially act by reducing the line tension Ao for the formation of free edges 
and seams. However, I do not see how this suggestion helps in understanding the asymmetry in the form of 
the scattering function I(q), or the apparent lack of divergence in 1(0). Concerning the last point the authors 
argue that the actual behavior of the membrane should be described by a more general model, in which edges 
and seams have different line tensions. However, the critical behavior of this model is still unknown. 

Plumber's nightmares 

I close by a very brief discussion of plumber's nightmare phases The free energy per unit volume of this 
phase as a function of the lattice constant s is given by 

t I \ ~ n T ° _l C 2 k B T - C 3 k(s) 

•F p „(s) ~ 61 1 ^ . (6.60) 

s s J 

The density of handles is C3/(47rs 3 ), and the constants Ci are of order unity. The lattice constant will 
presumably be of the order of the short-distance cutoff £q, at least if the "bare" gaussian rigidity kq is 
positive. In this case, by comparing this free energy to that of the droplet phase (which is zero on the scale 
t~o/£q) the transition to the plumber's nightmare phase will take place when C3K0 — Cik^T + C\Tq£q. We 
expect the transition from the droplet phase to be strongly first-order. On the other hand one can view the 
transition from the lamellar side to be analogous to a "crystallization" of the passage fluid considered, e.g., 
in ref. pS|]. In this case it should be weakly first order. 



Conclusions 

Unfortunately I have not been able to discuss in these notes some fascinating arguments, like the behavior 
of membranes with reduced symmetry (hexatic and cristalline) , or made of several components. But I had 
presumed too much on my speed in writing and I could not abuse of the editors' patience. I hope however 
that I have conveyed to you part of my fascination for this delicate field of Statistical Physics. 
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Figure 43: Axisymmetric vesicle shape: the z axis is the axis of rotational symmetry. 



A Differential equations for vesicle shapes 

The Euler-Lagrange equations for vesicle shapes cannot be analytically solved in general. We show in this 
Appendix how, limiting oneself to the consideration of axisymmetric shapes, one can transform them into a 
system of first-order ordinary differential equations, which can be solved numerically ((3^, |9^]). 
The geometry is defined in fig. (fl3]). Because of rotational symmetry, one has 



H = c m + c p , 



(A.l) 



where c m;P are the principal curvatures (the eigenvalues of the curvature matrix f2), corresponding respec- 
tively to the meridians (m) and the parallels (p). Denoting by ip the angle between the surface normal and 
the z-axis, we have 



sin-0(x) 



The angle ip(x) is related to the contour z{x) by 



it \ d ^' 
cosiplx)— — . 

ax 



dz 

— = — taxiipix). 
dx 

On the other hand, the contour z{x) and the volume V can be expressed in terms of x and c p : 

z(x) = z(0) - J dxxc p [l - (xcp) 2 ]^ 1 , 

dV = irx 3 c p [l - {xc p ) 2 ] 2 dx, 
dS = 2ttx [1 - {xc p ) 2 ]~ 2 ~ dx. 
It is useful to replace the independent variable x by the normalized surface area s using the relation 



dx = ±— [l — (xcp) 



21 2 



(A.2) 
(A.3) 

(A.4) 

(A.5) 
(A.6) 

(A.7) 



The variable s ranges from s = at the upper pole to a value s m at the dividing parallel of latitude, and 
then to s — 1 at the lower pole. One has 



da; J > 0, for < s < s m ; 
ds 1 < 0, for s m <s< 1. 



(A. 
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Introducing the notation / = x 2 the equations become: 



dc m 
ds 



dcp 
~d7 

ds 



= ± 



(! - / c p) 2 x { c p ( c p - c °) 



^ c 2 



±2 1 - ./(• 



2\5 u m 



±4 (l-/^) 5 



(A.9) 
(A.10) 
(A.ll) 



These equation have singularities for either f(s) — > or f(s)c 2 — > 1. At the poles of the vesicles we have 
/ = 0, but at the same time the difference c m — c p tends to zero. This yields the correct boundary condition 
at the poles: 



dcp 
ds 



poles 



1 dCm 

3 ds 



(A.I2) 



poles 



The singularity for f(s)c 2 — > 1 is approached when /(s) reaches an extremum, and should be removed by 



imposing that the factor in curly brackets in eq. (A. 10) vanishes. This allows one to join the ± branches of 
the equations by expanding in a Taylor series around the value s = s m where /(s) reaches the extremum. 

One can then obtain the three curves c m (s), c p (s), /(s), and the matching at those intermediate points s m 
where f(s)c 2 — ► 1 gives three nonlinear equations for the unknown boundary conditions c m (s=0), f(s=s m ), 
and c m (s—s m ). A further integration yields the actual vesicle shape z(x). 



B The Faddeev-Popov determinant 

Let us consider the integral 

1 = J d 2 xf(x), 

where the function f(x,y) is invariant upon rotations in the (x,y) plane: 

x — > xg — £1(9) ■ x, 

where 0,(9) is a rotation matrix: 



= 



cos f — sm t 
sin 9 cos 6 



(B.l) 



(B.2) 



(B.3) 



It is then convenient to integrate the function over one line, whose images upon rotation sweep the whole 
two-dimensional plane. This line can be represented by an equation in x: 



ip(x) = 0. 



(B.4) 



Given any vector x there will be a value 9* such that ip(x$*) — 0. Using the properties of the delta function, 
we obtain 



2tt 



d9S (ip(x e )) 



dip(x e ) 



09 



(B.5) 



Therefore 



I 



2 7T 



A9 I d 2 xS(ip(rg)) 



dip(x e ) 



09 



2ir / d 2 x6(?p(r)) 



dtj)(xe) 



09 



fif)- 



(B.6) 



The advantage of this expression is that all the quantities need only to be evaluated on (or near) the line 
ip(x) = 0. 
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The Faddeev-Popov determinant is the generalization of the factor introduced above to correct the effect 
of the integration over 9. We consider a general integral 



J dfj,{x) f(x), 



(B.7) 



where the measure d/j,(x) and the function f(x) are both invariant with respect to a TV-dimensional group 
G, whose elements g can be represented via its TV generators t a : 



9 = cxp 



We "fix the gauge", introducing N constraints ip a {x) = 0, i = 1, . . . , N. We then have 



N 



1= / d^(x)f(x) l[6(il> a (x)) |Det(J a/3 (a:))|volG, 



a=l 



where we have introduced the volume of the group G, and 

dlp a (x) 



Jat3{x) 



do? 



0=0 



In our case, the relevant group is generated by the infinitesimal diffeomorphisms of the form 

(7 — > o_ -|- e(a) . 
The corresponding transformations on the r read 

r — > f + e l dif. 

Consider now a deformation h with respect to a background configuration fb: 

= f (s.) + h(a). 

The normal gauge is realized by imposing the constraints 

i>i{s) = h(a) ■ dir (a) =0, i = 1,2. 

Let us compute how the constraints are affected by an infinitesimal diffeomorphism: 

6ipi(a) = S(r-r )- d l r = Sh ■ dif 
= e 3 djf-dif 
= e j (djf + djhj ■ difo 

= e J (gij + djh ■ <9V ) . 



(B.8) 



(B.9) 

(B.10) 

(B.11) 
(B-12) 
(B.I3) 
(B.14) 



(B.15) 



We see that this expression is purely local. Therefore, in this gauge, the Faddeev-Popov determinant reads 



A F p = Yl Dct 9ij is.) + dih(a) ■ djf {a) 



(B.I6) 



In the particular case of the Monge representation, dji{q) ■ djfo(a) = 0, so that the Faddeev-Popov deter- 
minant is trivial. 
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C One-loop calculation of the renormalization group flow 



In this Appendix I report the calculation of the renormalization group flow to one loop, following essentially 
plj , p8| , |5l| . We use the method of the background field. We take the reference configuration to be paramet- 
rically described by fo(a) and we sum over all configurations r(a) such that: (i) they are sufficiently "close" 
to the reference one (so that the deformation h(a) is small); (ii) the wavenumbers contained in the Fourer 
transform of h(a) belong to the infinitesimal shell A/s < \q\ < A. We write the effective Wilson hamiltonian 
H e s(s) as follows: 

F[r + h]\ 



Hefr(s) = — /cB^log / Vhexp 



dA\ ■ f , 



(C.l) 



where J is a reminder on the constraint on the fluctuation wavenumber, and A is the external field needed 

to set (r) x — 0. We can expand F[ro + h] in powers of h up to second order. The contribution of the first 
order will cancel with the second term of the rhs of eq. (C.l). We obtain 




Hes(s) 



where 5 2 F is the second variation of the Helfrich hamiltonian with respect to h. 
We now introduce the metric tensor g^ of the reference configuration r*o : 

9ij(s.) = dif (a) ■ djf (a). 



(C.2) 



(C.3) 



On the basis of gij we can define the covariant derivative Di . We can then express the curvature tensor fiy 
as follows: 



D l D j r . 



(C.4) 



One can check that fly is normal to the surface for any i, j. It can therefore represented by fJy 
where n is the normal unit vector. The Gaussian curvature K and the square of the mean curvature H 2 are 
then given by 



H 2 = fVlV; 



(C.5) 
(C.6) 



We use the normal gauge, and represent the deformation h by the scalar v, defined by h — vri. We then 
have: 



5 2 F 



AAv { k 



^(H 2 - AK) + •_'//<>"/),/), 



+t [-A + 2K]}p, 

where we have introduced the Laplacian (with respect to the metric tensor gij): 

A = DiD\ 



(C.7) 



(C.8) 



The gaussian curvature term drops out of the variation. We have neglected terms which are irrelevant in 
our approximations, being of higher order in the curvatures. When integrating upon the fluctuations, the 
Tr s log term in eq. (CI2) can be replaced by J daj s d 2 g(2-7r)~ 2 , and the operator A by q 2 . The result reads 



n cS (s) = F 



dA- 



x loe 



1 



k B TA 4 



2tt 

kA 4 + tA 2 



l K H 2 A 2 + —kK 
2 3 



(C.9) 
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We can read off this formula the effective parameters: 



«W - ^T^fe + i^)' (ai0) 



**« = K ~ £ 47 r l + (r/A') ' (C ' n) 



^ cff ( S ) = ^ + £ 5 ^ Br / (C.12) 
37r 1 + (r/ A^J 

However, this formula does not consider the contributions of the measure terms |fl3|j . Their effect is a shift 
in r off : 

r ^^ T eff = T ^- e ^ fc BTl0g('-^V (C.13) 

Taking into account this contribution and rescaling back the lengths we obtain the renormalized couplings, 
which obey the renormalization group equations 

dr _ A 



2r + - fcB Tlog(l + ^), (C.14) 
dn _ 3fc B T 1 

5s 4tt 1 + (t/A2)' 1 j 



5k _ 5fc B T 1 
S fo " 3tt 1 + (t/A*)' ' (,i 



D The Liouville model 

The Liouville model is defined by the hamiltonian 

W, gij] = 7o / d 2 aV5 + A J <ra x 7,g"i),r- 3 r. (D.l) 

The metric tensor g^ fluctuates in this model, and is not induced by the membrane configuration described 
by r. 

The Liouville model is best studied in the conformal gauge where the internal metric tensor g^ takes the 
form 

<?ii=e%. (D.2) 



It is always possible to construct locally such a coordinate system (see, e.g., [28, Sec. 2.8]). In this system 
the gaussian curvature reads 

1 / d 2 d ? 



K ^=-'r^[M + M)^- (D - 3) 

If we fix (f>, we can integrate over the Gaussian field r, obtaining an effective hamiltonian for the metric 
field (f>. This integration can be performed via the conformal anomaly [IOC, |28| ]. 

Consider a scalar field defined on a closed surface S, and associate to it a hamiltonian of the form 

n [if>] = \J dV^'ftH^- (d.4) 

Integrating by parts we obtain 

Ho bP] = \J dV VflW- AM (D.5) 

where A is the scalar Laplacian in the metric gij: A = g v DiDj. Di is the covariant derivative with respect 
to a 1 . In the conformal gauge, the Laplacian is simply given by 
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Integrating over the Gaussian field ip yields formally 

Z = j Vip e- Ho [if;} oc [Det(-A)] _ * . 
One can express it formally in terms of the eigenvalues {— A&} of the Laplacian: 



Det(-A) = J] A fe = exp j £ log X k J 

fe V k ) 



(D.7) 



(D. 



However, this expression is ill-defined. On the one hand, the constant function belongs to the eigenvalue zero 
(yielding an infinite contribution). On the other hand, the sum diverges at large k. One has to introduce an 
ultraviolet (large wavenumber) regulator A ~ 7r/aoi to suppress this divergence. This regulator should be 
introduced in a reparametrization-invariant way, which is obtained via the heat-kernel regularization. 
The "heat kernel" G(a,a';t) is the solution of the differential equation 



^G( 2)2 , ;i)=A CT G( £ , 2 ';t), 



satisfying the initial condition 



G(a,2';i=0) 



5{g_-a!). 



(D.9) 



(D.10) 



The extra factor on the rhs is introduced to maintain reparametrization invariance. At short times, G(a_, a'; t) 
looks locally like a Gaussian, with a width proportional to ^ft, like a diffusion kernel. However, this diffusion 
takes place in a curved space. It will therefore be slower in a space of positive curvature and faster in a space 
of negative curvature. The asymptotic expansion of G(<7, a; t) for t — ► + is indeed given by [ pl| 

1 Tl K(a) 
4^ t + 3 



G(a,a; t) 



0{t) 



(D.ll) 



We can now define the regularized determinant of the Laplacian as follows: 



logDet'(-Ae) = Tr'log(-A £ 



Tr' 



^Tr' (e* A ) 



(D.12) 



The parameter e ~ A 2 acts as a cutoff. The prime indicates that the zero eigenvalue has been suppressed. 
By using the expression ( D.10| ) we obtain 



Tr'log(-A e ) ~ --!- / d 2 a^j- 1 ^- [ & 2 JgK + finite terms. 
4ttc Js 12n J s 



(D.13) 



Note that the coefficient of 1/e depends on the cutoff procedure used, while the coefficient of loge does not. 
The finite contribution can be computed in the following way. Consider a local conformal rescaling 



^■fe)^i(£)e^=<4-. 
Under such a transformation, the Laplacian changes as 



(D.14) 
(D.15) 



Therefore 



S(p (a) 



Tr'log(-A e ) = 



00 dt 
— Tr 



— e 

Sip 



dt 



/°° A 



1 K(a) 
47re 12n 



O(e) 



(D.16) 
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This is the conformal anomaly | IOC ] . This result can also be written in the remarkable form: 



Sip{a 



■Tr'log(-A e ) = 



VgW)K(a)G(a, a') y/g^K^) 



(D.17) 



where G(q_, a') is the inverse of the Laplacian in the metric gij. Therefore one can integrate out the variation 
in v?(<l)j starting, e.g., from a constant metric conformally equivalent to gij, and obtain 



Trlog(-A e (. 9 )) = -L f d 2 a^/gjaj 

** s 



-T[g°] + 0(e), 



(D.18) 



where ^[ff ] depends only on the conformal class of the metric g, and contains the logarithmically divergent 
part obtained in (D.13). In the conformal gauge, this expression becomes simply 



ille 



48tt 



dV 



12 



(D.19) 



The finite part in eq. (D.1S) is called the Liouville action [100]. The Liouville model (D.l) involves d scalar 
fields, corresponding to the components of the field f. On the other hand, the Faddeev-Popov determinant 
(App. B) is nontrivial in the conformal gauge. Its contribution turns out to be also proportional to the 
Liouville action, but with a different coefficient. The effective hamiltonian for d> reads 



48tt J s 



(D.20) 



where fj, 2 is a "mass" which represents the effect of the surface tension. 



References 

[1] D. Andclman, Electrostatic properties of membranes: The Poisson-Boltzmann theory, in R. Lipowsky, 
E. Sackmann (eds.) Membranes: Their structure and conformation (Elsevier, 1994). 



[2 
[3 

[4: 

[5 

[6 

[7: 

[8 
[9 

[io; 
[ii 

[12 



J. R. Banavar, A. Maritan, A. Stella, Phys. Rev. A43 5752 (1991). 
A. Baumgartner, Physica A190 63 (1992). 

K. Berndl, J. Kas, R. Lipowsky, E. Sackmann, U. Seifert, Europhys. Lett. 13 659 (1990). 
M. Bloom, E. Evans, O. G. Mouritsen, Quarterly Review of Biophysics 24 293 (1991). 
M. Blume, V. J. Emery, R. B. Griffiths, Phys. Rev. 67 3194 (1991). 
D. H. Boal, M. Rao, Phys. Rev. A 46 3037 (1992). 

D. H. Boal, U. Seifert, A. Zilker, Phys. Rev. Lett. 69 3405 (1992). 
K. A. Brakke, J. Exp. Math. 1 141 (1992). 

F. Brochard, P.-G. de Gennes, P. Pfeuty, J. Phys. France 37 1099 (1976). 
F. Brochard, J. F. Lennon, J. Phys. France 36 1035 (1975). 

E. Browicz, Zbl. Med. Wiss. 28 625 (1890). 



62 



[13] W. Cai, T. C. Lubensky, P. Nelson, T. Powers, J. Phys. II France 4 931 (1994). 
[14] P. B. Canham, J. theor. Biol. 26 61 (1970). 

[15] L. Cantii, M. Corti, M. Musolino, P. Salina, Europhys. Lett. 13 561 (1991). 
[16] L. Cantii, M. Corti, E. Del Favcro, A. Raudino, J. Phys. II France 4 1585 (1994). 
[17] See the lectures by J. Cardy in this School. 
[18] M. E. Cates, Europhys. Lett 8 719 (1988). 

[19] M. E. Cates, D. Roux, D. Andelman, S. Milner, S. Safran, Europhys. Lett. 5 733 (1988); Erratum ibid. 
7 94 (1988). 

[20] M. E. Cates, D. Roux, J. Phys. Cond. Matter 2 SA339 (1990). 
[21] G. Cevc (ed.), Phospholipids handbook (New York: Dekker, 1993). 

[22] G. Cevc, D. Marsh, Phospholipid bilayers: Physical principles and models (New York: Wiley, 1987). 

[23] D. Chapman, et al., Chem. Phys. Lipids 1 445 (1967). 

[24] C. Coulon, D. Roux, A. M. Bellocq, Phys. Rev. Lett. 66 1709 (1991). 

[25] C. Coulon, D. Roux, M. E. Cates, Phys. Rev. Lett. 67 3194 (1991). 

[26] B. Dammann, H. C. Fogedby, J. H. Ipsen, C. Jeppescn, Preprint Aarhus IFA-94/11. 

[27] F. David, Europhys. Lett. 2 577 (1986). 

[28] F. David, in: D. Nelson, T. Piran, S. Weinberg (eds.), Statistical mechanics of membranes and surfaces 
(Singapore: World Scientific, 1989) p. 157ff. 

[29] F. David, J. Phys. France Suppl. n. 23 51 C7-115 (1990). 

[30] F. David, E. Guitter, Nucl. Phys. B295 332 (1988). 

[31] F. David, S. Leibler, J. Phys. II France 1 959 (1991). 

[32] T. Davis, et al., in: J. Meunier, D. Langevin, N. Boccara (eds.), Physics of amphiphilic layers (Berlin: 
Springer, 1987) p. 310. 

[33] H. J. Dculing, W. Helfrich, Bwphys. J. 16 861 (1976). 

[34] H. J. Deuling, W. Helfrich, J. Phys. France 37 1335 (1976). 

[35] J.-M. Drouffc, J.-B. Zubcr, Phys. Rep. 102 1 (1983). 

[36] B. Durhuus, J. Frohlich, T. Jonsson, Nucl. Phys. B225 185 (1983). 

[37] H. P. Duwc, J. Kas, E. Sackmann, J. Phys. France 51 945 (1990). 

[38] A. Elgsaeter, B. T. Stokke, A. Mikkelsen, D. Branton, Science 234 1218 (1986). 

[39] E. Evans, W. Rawicz, Phys. Rev. Lett. 64 2094 (1990). 

[40] M. Filali, G. Porte, J. Appell, P. Pfeuty, J. Phys. II France 4 349 (1994). 

[41] M. E. Fisher, Phys. Rev. 176 257 (1968). 

[42] D. Forster, Phys. Lett. 114A 115 (1986). 

[43] P. Fromherz, C. Rocker, D. Ruppcr, Faraday Discuss. Chem. Soc. 81 39 (1986). 



63 



[44] D. Gazeau, A. M. Bcllocq, D. Roux, T. Zcmb, Europhys. Lett. 9 447 (1989). 

[45] P.-G. dc Gcnncs, C. Taupin, J. Phys. Chem. 86 2294 (1982). 

[46] R. E. Goldstein, S. Lcibler, Phys. Rev. Lett. 61 2213 (1988). 

[47] R. E. Goldstein, S. Lcibler, Phys. Rev. A40 1025 (1989). 

[48] L. Golubovic, Phys. Rev. E50 R2419 (1994). 

[49] G. Gompper, D. M. Kroll, Europhys. Lett. 9 59 (1989) 

[50] G. Gompper, D. M. Kroll, Europhys. Lett. 19 581 (1992); Phys. Rev. A46 7466 (1992). 
[51] P. Greiner, Arch. Ration. Mech. and Analysis 41 163 (1971). 

[52] E. Guittcr, Physique statistique et proprietes critiques des membranes, These Paris VI (1989). 

[53] W. Harbich, W. Hclfrich, Z. Naturforsch. 34a 1063 (1979). 

[54] W. Harbich, R. M. Servuss, W. Hclfrich, Z. Naturforsch. 33a 1013 (1978). 

[55] W. Hclfrich, Z. Naturforsch. 28c 693 (1973). 

[56] W. Hclfrich, Z. Naturforsch. 33a 305 (1978). 

[57] W. Hclfrich, Amphiphilic mesophases made of defects, in: R. Balian et al. (eds.), Physics of defects, 
Les Houches XXXV (Amsterdam: North-Holland, 1981). 

[58] W. Hclfrich, J. Phys. France 46 1263 (1985). 

[59] W. Helfrich, J. Phys. France 47 321 (1986). 

[60] W. Helfrich, J. Phys. LL France 3 385 (1993). 

[61] P. Hcrvc, D. Roux, A.-M. Bcllocq, F. Nallet, T. Gulik-Krzywicki, J. Phys. II France 3 1255 (1993). 
[62] H. Hoffmann, C. Thunig, M. Valiente, Coll. Surf. Sci. 67 223 (1992). 

[63] H. Hoffmann, C. Thunig, U. Munkcrt, H. W. Meyer, W. Richtcr, Langmuir 8 2629 (1992). 

[64] L. Hsu, R. Kusner, J. Sullivan, J. Exp. Math. 1 191 (1992). 

[65] D. Huse, S. Leibler, J. Phys. France 49 605 (1988). 

[66] D. A. Huse, S. Leibler, Phys. Rev. Lett. 66 437 (1991). 

[67] M. J. Janiak, D. M. Small, G. G. Shipley, J. Biol. Chem. 254 6068 (1979). 

[68] G. A. Jongeward, J. D. Stack, C. Jayaprakash, Phys. Rev. D21 3360 (1980). 

[69] F. Jiilicher, U. Seifert, R. Lipowsky, J. Phys. II France 3 1681 (1993). 

[70] J. Kas, E. Sackmann, Biophys. J. 60 825 (1991). 

[71] J. Kas, E. Sackmann, R. Podgornik, S. Svetina, B. Zeks, J. Phys. II France 3 631 (1993). 
[72] C. L. Khetrapal, et al., Lyotropic liquid crystals (Berlin: Springer, 1975). 
[73] H. Klcincrt, Phys. Lett. A114 263 (1986). 

[74] R. D. Romberg, H. M. McConncll, Biochemistry 10 1111 (1971). 

[75] L. D. Landau, E. M. Lifshitz, Theory of elasticity (Oxford: Pcrgamon, 1980). 

[76] F. Larche, J. Appcll, G. Porte, P. Bassereau, J. Marignan, Phys. Rev. Lett. 56 1700 (1986). 



64 



[77] S. Lcibler, J. Phys. France 47 507 (1986). 

[78] S. Leibler, in: D. Nelson, T. Piran, S. Weinberg (eds.), Statistical mechanics of membranes and surfaces 
(Singapore: World Scientific, 1989) p. 45ff. 

[79] S. Leibler, R. R. P. Singh, M. E. Fisher, Phys. Rev. Lett. 59 1989 (1987). 

[80] R. Lipowsky, Nature 349 475 (1991). 

[81] R. Lipowsky, S. Lcibler, Phys. Rev. Lett. 56 2541 (1986); erratum, ibid. 56 2541 (1986). 

[82] R. Lipowsky, D. Richtcr, K. Kremer (eds.), The structure and conformation of amphiphilic membranes 
(Berlin: Springer, 1992). 

[83] N. D. Mcrmin, H. Wagner, Phys. Rev. Lett. 17 1133 (1966). 

[84] J. Meunier, D. Langevin, N. Boccara (eds), Physics of amphiphilic layers (Berlin: Springer, 1987). 
[85] L. Miao, U. Seifert, M. Wortis, H.-G. Dobcrcincr, Phys. Rev. E49 5389 (1994). 

[86] X. Michalct, Etude experimental de vesicules phospholipidiques de genre topologique non spherique, 
These Paris VII (1994). 

[87] X. Michalet, D. Bcnsimon, B. Fourcade, Phys. Rev. Lett. 72 168 (1994). 

[88] S. T. Milncr, D. Roux, J. Phys. I France 2 1741 (1992). 

[89] D. C. Morse, Phys. Rev. E50 R2423 (1994). 

[90] D. C. Morse, S. T. Milncr, Europhys. Lett. 26 565 (1994). 

[91] M. Mutz, W. Helfrich, J. Phys. France 51 991 (1990). 

[92] M. Mutz, D. Bcnsimon, Phys. Rev. A43 4525 (1991). 

[93] D. Nelson, T. Piran, S. Weinberg (eds.), Statistical mechanics of membranes and surfaces (Singapore: 
World Scientific, 1989). 

[94] R. H. Pearson, I. Pascher, Nature 281 499 (1979). 

[95] L. Pcliti, S. Leibler, Phys. Rev. Lett. 54 1690 (1985). 

[96] M. A. Peterson, J. Appl. Phys. 57 1739 (1985). 

[97] M. A. Peterson, J. Math. Phys. 26 711 (1985). 

[98] M. A. Peterson, H. Strcy, E. Sackmann, J. Phys. II France 2 1273 (1992). 
[99] A. G. Petrov, J. Bivas, Progress Surf. Sci. 18 389 (1984). 
[100] A. M. Polyakov, Phys. Lett. B103 207 (1981). 

[101] G. Porte, J. Appell, P. Bassereau, J. Marignan, J. Phys. France 49 511 (1988). 

[102] G. Porte, J. Marignan, P. Bassereau, R. May, J. Phys. France 49 511 (1988). 

[103] T. Powers, Seminar on "Measure factors for fluctuating membranes" , this School. 

[104] D. Roux, M. E. Cates, U. Olsson, R. G. Ball, F. Nallet, A. M. Bellocq, Europhys. Lett. 10 229 (1990). 

[105] D. Roux, C. Coulon, M. E. Cates, J. Phys. Chem. 96 4174 (1992). 

[106] E. Sackmann ct al., in Liquid crystals of one- and two-dimensional order (Berlin: Springer, 1980). 
[107] U. Seifert, J. Phys. A: Math. Gen. 24 L573 (1991). 



65 



[108] U. Seifert, K. Berndl, R. Lipowsky, Phys. Rev. A 42 4768 (1991). 

[109] B. L. Silver, The physical chemistry of membranes (Boston: Allen & Unwin, 1985). 

[110] B. D. Simons, M. E. Cates, J. Phys. II France 2 1439 (1992). 

[Ill] S. J. Singer, Science 175 720 (1972). 

[112] G. S. Smith, E. B. Sirota, C. R. Safinya, N. A. Clark, Phys. Rev. Lett. 60 813 (1988). 

[113] R. Strey, W. Jahn, G. Porte, P. Bassereau, Langmuir 6 1635 (1990). 

[114] L. Strycr, Biochemistry 3rd ed. (New York: Freeman, 1988). 

[115] S. Svetina, A. Ottova-Lietmannova, R. Glascr, J. theor. Biol. 94 13 (1982). 

[116] I. Szlcifcr, D. Kramer, A. Ben-Shaul, D. Roux, W. M. Gelbart, Phys. Rev. Lett. 60 1966 (1988). 
[117] Y. Talmon, S. Pragcr, J. Chem. Phys. 69 2984 (1978). 
[118] D. C. Wack, W. W. Webb, Phys. Rev. Lett. 61 1210 (1988). 

[119] T. J. Willmorc, Total curvature in Riemannian geometry (New York: Horwood, 1982). 



